(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 

(19) World Intellectual Property Organization 

International Bureau 

(43) International Publication Date r>r"T ^ Interna ti° nal Publication Number 

10 May 2007 (10.05.2007) WO 2007/053486 A2 





(51) International Patent Classification: 
G06F 79/00(2006.01) 

(21) International Application Number: 

PCT/US2006/042140 

(22) International Filing Date: 30 October 2006 (30.10.2006) 



(25) Filing Language: 

(26) Publication Language: 



(30) Priority Data: 

60/730,882 
60/732,154 
60/737,744 
6(V758,528 
60/780,518 
60/788,694 
60/812,590 
60/815,253 



28 October 2005 
2 November 2005 
18 November 2005 
13 January 2006 
9 March 2006 
4 April 2006 
12 June 2006 
21 June 2006 



English 
English 



(28.10.2005) US 

(02.11.2005) US 

(18.11.2005) US 

(13.01.2006) US 
(09.03.2006) US 
(04.04.2006) US 
(12.06.2006) US 
(21.06.2006) US 



(71) Applicant (for all designated States except US): BLACK- 
LIGHT POWER, INC. [US/US] ; 493 Old Trenton Road, 
Cranbury,NJ 08512 (US). 

(72) Inventor; and 

(75) Inventor/Applicant (for US only): MILLS, Randeli, L. 
[US/US]; 27 Rosedale Road, Princeton, NJ 08540 (US). 



(74) Agent: MELCHER, Jeffrey, S.; Manelli Denison & Sel- 
ler, PLLC, 2000 M Street, N.W., 7th Floor, Washington, 
DC 20036-3307 (US). 

(81) Designated States ( unless otherwise indicated, for every 
kind of national protection available): AE, AG, AL, AM, 
AT, AU, AZ, BA, BB, BG, BR, B W, BY, BZ, CA, CH, CN, 
CO, CR, CU, CZ, DE, DK, DM, DZ, EC, EE, EG, ES, FI, 
GB, GD, GE, GH, GM, GT, HN, HR, HU, ID, IL, IN, IS, 
JP, KE, KG, KM, KN, KP, KR, KZ, LA, LC, LK, LR, LS, 
LT, LU, LV, LY, MA, MD, MG, MK, MN, MW, MX, MY, 
MZ, NA, NG, NT, NO, NZ, OM, PG, PH, PL, PT, RO, RS, 
RU, SC, SD, SE, SG, SK, SL, SM, SV, SY, TJ, TM, TN, 
TR, TT, TZ, UA, UG, US, UZ, VC, VN, ZA, ZM, ZW. 

(84) Designated States (unless otherwise indicate a\ for every 
kind of regional protection available): ARIPO (BW, GH, 
GM, KE, LS, MW, MZ, NA, SD, SL, SZ, TZ, UG, ZM, 
ZW), Eurasian (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), 
European (AT, BE, BG, CH, CY, CZ, DE, DK, EE, ES, FI, 
FR, GB, GR, HU, IE, IS, IT, LT, LU, LV, MC, NL, PL, PT, 
RO, SE, SI, SK, TR), OAPI (BF, BJ, CF, CG, CI, CM, GA, 
GN, GQ, GW, ML, MR, NE, SN, TD, TG). 

Published: 

— without international search report and to be republished 
upon receipt of that report 

For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCT Gazette. 



(54) Title: SYSTEM AND METHOD OF COMPUTING AND RENDERING THE NATURE OF POLYATOMIC MOLECULES 
AND POLYATOMIC MOLECULAR IONS 



Start 



Functional group data files. 
Molecules data files 



Process input files according to 
their format 



Make functional -group objects from 
Functional groups data Ales 



i SMILES input, break it down 
into fiff)Ctior^l £Foups 

i 1 1 ■ " 



Calculate geometric parameters of 
Molecules depicted in molecules 
data files from functional- group 

objects, end make molecule objects 





< 
00 

m 
m 

(57) Abstract: A method and system of physically solving the charge, mass, and current density functions of polyatomic molecules, 
polyatomic molecular ions, diatomic molecules, molecular radicals, molecular ions, or any portion of these species using Maxwell's 
equations and computing and rendering the physical nature of the chemical bond using the solutions. The results can be displayed 

^ on visual or graphical media. The display can be static or dynamic such that electron motion and specie's vibrational, rotational, and 
translational motion can be displayed in an embodiment. The displayed information is useful to anticipate reactivity and physical 

^ properties. The insight into the nature of the chemical bond of at least one specie can permit the solution and display of those of 
other species to provide utility to anticipate their reactivity and physical properties. 
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SYSTEM AND METHOD OF COMPUTING AND RENDERING THE NATURE OF 
POLYATOMIC MOLECULES AND POLYATOMIC MOLECULAR IONS 

This application claims priority to U.S. Application Nos.: 60/730,882, filed October 28, 
2005; 60/732,154, filed November 2, 2005; 60/737,744, filed November 18, 2005; 
60/758,528, filed January 13, 2006; 60/780,518, filed March 9, 2006; 60/788,694, filed April 
4, 2006; 60/812,590, filed June 12, 2006; and 60/815,253, June 21, 2006, the complete 
disclosures of which are incorporated herein by reference. 

Field of the Invention: 

This invention relates to a system and method of physically solving the charge, mass, 

and current density functions of polyatomic molecules, polyatomic molecular ionsj diatomic 

* • - ' * 

' : « * . , * • ' * ■ ■ - ' • 

molecules, molecular radicals, molecular ions, or any portion of these species, and computing 
and rendering the nature of these species using the solutions. The results can be displayed on 
visual or graphical media. The displayed information provides insight into the nature of these 
species and is useful to anticipate their reactivity, physical properties, and spectral absorption 
and emission, and permits the solution and display of other species. 

Rather than using postulated unverifiable theories that treat atomic particles as if they 
were not real, physical laws are now applied to atoms and ions. In an attempt to provide 
some physical insight into atomic problems and starting with the same essential physics as 
Bohr of the e~ moving in the Coulombic field of the proton with a true wave equation, as 
opposed to the diffusion equation of Schitfdinger, a classical approach is explored which 
yields a model that is remarkably accurate and provides insight into physics on the atomic 
level. The proverbial view deeply seated in the wave-particle duality notion that there is no 
large-scale physical counterpart to the nature of the electron is shown not to be correct 
Physical laws and intuition may be restored when dealing with the wave equation and 
quantum atomic problems. 

Specifically, a theory of classical quantum mechanics (CQM) was derived from first 
principles as reported previously [reference Nos. 1-8] that successfully applies physical laws 
to the solution of atomic problems that has its basis in a breakthrough in the understanding of 
the stability of the bound electron to radiation. Rather than using the postulated Schrtfdinger 
boundary condition: " *F -> 0 as r -» oo which leads to a purely mathematical model of the 
electron, the constraint is based on experimental observation. Using Maxwell's equations* the 
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classical wave equation is solved with the constraint that the bound n = 1 -state electron 
cannot radiate energy. Although it is well known that an accelerated point particle radiates, 
an extended distribution modeled as a superposition of accelerating charges does not have to 
radiate. A simple invariant physical model arises naturally wherein the predicted results are 
extremely straightforward and internally consistent requiring minimal math, as in the case of 
the most famous equations of Newton, Maxwell, Einstein, de Broglie, and Planck on which 
the model is based. No new physics is needed; only the known physical laws based on direct 
observation are used. 

Applicant's previously filed WO2005/067678 discloses a method and system of 
physically solving the charge, mass, and current density functions of atoms and atomic ions 
and computing and rendering the nature of these species using the solutions. . The complete 
disclosure of this published PCT application is incorporated herein by reference. 

Applicant's previously filed WO2005/1 16630 discloses a method and system of 
physically solving the charge, mass, and current density functions of excited states of atoms 
and atomic ions and computing and rendering the nature of these species using the solutions. 
The complete disclosure of this published PCT application is incorporated herein by 
reference. 

Applicant's previously filed U.S. Published Patent Application No. 20050209788A1, 
relates to a method and system of physically solving the charge, mass, and current density 
functions of hydrogen-type molecules and molecular ions and computing and rendering the 
nature of the chemical bond using the solutions. The complete disclosure of this published 
application is incorporated herein by reference. 

Background of the Invention 

The old view that the electron is a zero or one-dimensional point in an all-space 
probability wave function is not taken for granted. The theory of classical quantum 
mechanics (CQM), derived from first principles, must successfully and consistently apply 
physical laws on all scales [1-8]. Stability to radiation was ignored by all past atomic models. 
Historically, the point at which QM broke with classical laws can be traced to the issue of 
nonradiation of the one electron atom. Bohr just postulated orbits stable to radiation with the 
further postulate that the bound electron of the hydrogen atom does not obey Maxwell's 
equations — rather it obeys different physics [1-12]. Later physics was replaced by "pure 
mathematics" based on the notion of the inexplicable wave-particle duality nature of electrons 
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which lead to the Schrtidinger equation wherein the consequences of radiation predicted by 
Maxwell's equations were ignored. Ironically, Bohr, Schrddinger, and Dirac used the 
Coulomb potential, and Dirac used the vector potential of Maxwell's equations. But, all 
ignored electrodynamics and the corresponding radiative consequences. Dirac originally 
attempted to solve the bound electron physically with stability with respect to radiation 
according to Maxwell's equations with the further constraints that it was relativistically 
invariant and gave rise to electron spin [13]. He and many founders of QM such as 
Sommerfeld, Bohm, and Wein stein wrongly pursued a planetary model, were unsuccessful, 
and resorted to the current mathematical-probability-wave model that has many problems [9- 
16]. Consequently, Feynman for example, attempted to use first principles including 
Maxwell's equations to discover new physics, to replace quantum mechanics [17]. 

Physical laws may indeed be the root of the observations thought to be "purely 
quantum mechanical", and it was a mistake to make the assumption that Maxwell's 
electrodynamic equations must be rejected at the atomic level. Thus, in the present approach, 
the classical wave equation is solved with the constraint that a bound n = 1 -state electron 
cannot radiate energy. 

Herein, derivations consider the electrodynamic effects of moving charges as well as 
the Coulomb potential, and the search is for a solution representative of the electron wherein 
there is acceleration of charge motion without radiation. The mathematical formulation for 
zero radiation based on Maxwell's equations follows from a derivation by Haus [18], The 
function that describes the motion of the electron must not possess spacetime Fourier 
components that are synchronous with waves traveling at the speed of light. Similarly, 
nonradiation is demonstrated based on the electron's electromagnetic fields and the Poynting 
power vector. 

It was shown previously [1-8] that CQM gives closed form solutions for the atom 
including the stability of the n - 1 state and the instability of the excited states, the equation 
of the photon and electron in excited states, and the equation of the free electron and photon, 
which predict the wave particle duality behavior of particles and light. The current and 
charge density functions of the electron may be directly physically interpreted. For example, 
spin angular momentum results from the motion of negatively charged mass moving 
systematically, and the equation for angular momentum, r x p, can be applied directly to the 
wave function (a current density function) that describes the electron. The magnetic moment 
of a Bohr magneton, Stern Gerlach experiment, g factor, Lamb shift, resonant line width and 
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shape, selection rules, correspondence principle, wave particle duality, excited states, reduced 
mass, rotational energies, and momenta, orbital and spin splitting, spin-orbital coupling, 
Knight shift, and spin-nuclear coupling, and elastic electron scattering from helium atoms, are 
derived in closed-form equations based on Maxwell's equations. The calculations agree with 
experimental observations. 

The Schrddinger equation gives a vague and fluid model of the electron. Schrodinger 
interpreted e*¥ * (xy¥(x) as the charge-density or the amount of charge between x and 
x + dx (T * is the complex conjugate of *F ). Presumably, then, he pictured the electron to 

■ 

be spread over large regions of space. After SchrSdinger's interpretation, Max Born, who 
was working with scattering theory, found that this interpretation led to inconsistencies, and 
he replaced the Schrfldinger interpretation with the probability of finding the electron . 
between x end.x + dx as 



Bora's inteipretation is generally accepted. Nonetheless, interpretation of the wave function 
is a never-ending source of confusion and conflict Many scientists have solved this problem 
by conveniently adopting the SchrSdinger interpretation for some problems and the Born 
interpretation for others. This duality allows the electron to be everywhere at one time— yet 
have no volume. Alternatively, the electron can be viewed as a discrete particle that moves 
here and there (from r = 0 to r = oo ), and * gives the time average of this motion. 

In contrast to the failure of the Bohr theory and the nonphysical, adjustable-parameter 

i 

approach of quantum mechanics, multielectron atoms [1,4] and the nature of the chemical 
bond [1, 5] are given by exact closed-form solutions containing fundamental constants only. 
Using the nonradiative wave equation solutions that describe the bound electron having 
conserved momentum and energy, the radii are determined from the force balance of the 
electric, magnetic, and centrifugal forces that corresponds to the minimum of energy of the 
system. The ionization energies are then given by the electric and magnetic energies at these 
radii. The spreadsheets to calculate the energies from exact solutions of one through twenty- 
electron atoms are given in '06 Mills GUT [1] and are available from the internet [19]. For 
400 atoms and ions, as well as hundreds of molecules, the agreement between the predicted 
and experimental results is remarkable. 

The background theory of classical quantum mechanics (CQM) for the physical 
solutions of atoms and atomic ions is disclosed in R. Mills, The Grand Unified Theory of 
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Classical Quantum Mechanics, January 2000 Edition, BlackLight Power, Inc., Cranbury, 
New Jersey, (" '00 Mills GUT ,! )> provided by BlackLight Power, Inc., 493 Old Trenton Road, 
Cranbury, NJ, 08512; R. Mills, The Grand Unified Theory of Classical Quantum Mechanics, 
September 2001 Edition, BlackLight Power, Inc., Cranbury, New Jersey, Distributed by 
Amazon.com (" '01 Mills GUT), provided by BlackLight Power, Inc., 493 Old Trenton 
Road, Cranbury, NJ, 08512; R. Mills, The Grand Unified Theory of Classical Quantum 
Mechanics, July 2004 Edition, BlackLight Power, Inc., Cranbury, New Jersey, (" '04 Mills 
GUT"), provided by BlackLight Power, Inc., 493 Old Trenton Road, Cranbury, NJ, 08512; R. 
Mills, The Grand Unified Theory of Classical Quantum Mechanics, January 2005 Edition, 
BlackLight Power, Inc., Cranbury, New Jersey, (" f 05 Mills GUT), provided by BlackLight 
Power, Inc., 493 Old Trenton Road, Cranbury, N J, 085 1 2 (posted at 
www.blacklightrower.comk R. L. Mills, "The Grand Unified Theory of Classical Quantum 
Mechanics", June 2006 Edition, Cadmus Professional Communications-Science Press 
Division, Ephrata, PA, ISBN 0963517171, Library of Congress Control Number 
2005936834, (" '06 Mills GUT"), provided by BlackLight Power, Inc., 493 Old Trenton 
Road, Cranbury, NJ, 085 1 2 (posted at 

http://www.blacklightpower.com/bookdownload.shtnil) : in prior PCT applications 
PCT/US02/35872; PCT/US02/06945; PCT/US02/06955; PCT/US01/09055; PCT/US01/ 
25954; PCT/US00/20820; PCT/US00/20819; PCT/US00/09055; PCT/US99/17171; 
PCTAJS99/17129; PCT/US 98/22822; PCT/US98/14029; PCT/US96/07949; 
PCT/US94/02219; PCT/US9 1/08496; PCT/US90/01998; and PCT/US89/05037 and U.S. 
Patent No. 6,024,935; the entire disclosures of which are all incorporated herein by reference 
(hereinafter "Mills Prior Publications"). 

The following list of references, which are also incorporated herein by reference in 
their entirety, are referred to in the above sections using [brackets]: 

1 . R. L. Mills, "The Grand Unified Theory of Classical Quantum Mechanics", June 2006 
Edition, Cadmus Professional Communications-Science Press Division, Ephrata, PA, 
ISBN 0963517171, Library of Congress Control Number 2005936834; posted at 
http://www.blackliditoower.com/bookdowiiload.shto 

2. R. L. Mills, "Classical Quantum Mechanics", Physics Essays, Vol. 1 6, No. 4, December, 
(2003), pp. 433-498; posted with spreadsheets at 

www.blacklightpo wer . com/techpapers . shtml. 
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3. R. Mills, "Physical Solutions of the Nature of the Atom, Photon, and TTieir Interactions to 
Form Excited and Predicted Hydrino States", submitted. 

4. R. L. Mils, "Exact Classical Quantum Mechanical Solutions for One- Through Twenty- 
Electron Atoms", in press, posted with spreadsheets at 
http://www.blacklightpower.com/techpapers.shtml. 

5. R. L. Mills, "The Nature of the Chemical Bond Revisited and an Alternative Maxwellian 
Approach", Physics Essays, Vol. 17, (2004), pp. 342-389, posted with spreadsheets at 
http://www.blacklightpower.com/techpapers.shtml. 

6. R. L. Mills, "Maxwell's Equations and QED: Which is Fact and Which is Fiction", in 
press, posted with spreadsheets at http://www.blacklightpower.com/techpapers.ditml. 

7. R. L. Mills, "Exact Classical Quantum Mechanical Solution for Atomic Helium Which 
Predicts Conjugate Parameters from a Unique Solution for the First Time", submitted, 
posted with spreadsheets at http://www.blackhghlpower,com/theory/theory.shtml.. 

8. R. Mills, "The Grand Unified Theory of Classical Quantum Mechanics", Int. J. Hydrogen 
Energy, Vol. 27, No. 5, (2002), pp. 565-590. 

9. R. L. Mills, "The Fallacy of Feynman's Argument on the Stability of the Hydrogen Atom 
According to Quantum Mechanics", Annales de la Fondation Louis de Broglie, Vol. 30, 
No. 2, (2005), pp. 129-151, posted at http://www.blacklightpower.com/techpapers.shtml. 

10. R. Mills, The Nature of Free Electrons in Superfluid Helium—a Test of Quantum 
Mechanics and a Basis to Review its Foundations and Make a Comparison to Classical 
Theory, Int. J. Hydrogen Energy, Vol. 26, No. 10, (2001), pp. 1059-1096. 

11. R- Mills, "The Hydrogen Atom Revisited", Int. J. of Hydrogen Energy, Vol. 25, Issue 12, 
December, (2000), pp. 1171-11 83. 

12. F. Laloe, Do we really understand quantum mechanics? Strange correlations, paradoxes, 
and theorems, Am. J. Phys. 69 (6), June 2001, 655-701 . 

13. P. Pearle, Foundations of Physics, "Absence of radiationless motions of relativistically 
rigid classical electron", Vol. 7, Nos. 11/12, (1977), pp. 931-945. 

14. V. F. Weisskopf, Reviews of Modern Physics, Vol. 21, No. 2, (1949), pp. 305-315. 

15. H. Wergeland, "The Klein Paradox Revisited", Old and New Questions in Physics, 
Cosmology, Philosophy, and Theoretical Biology, A. van der Merwe, Editor, Plenum 
Press, New York, (1983), pp. 503-515. 

16. A. Einstein, B. Podolsky, N. Rosen, Phys. Rev., Vol. 47, (1935), p. 777. 

17. F. Dyson, "Feynman's proof of Maxwell equations", Am. J. Phys., Vol. 58, (1990), pp. 
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209-21 L 

18. Haus, H. A., "On the radiation from point charges", American Journal of Physics, 54, 
(1986), pp. 1126-1129. 

19. http://www.blacklightpowercom/new.shtml. 

SUMMARY OF THE INVENTION 

The present invention, an exemplary embodiment of which is also referred to as 
Millsian software, stems from a new fundamental insight into the nature of the atom. 
Applicant's new theory of Classical Quantum Mechanics (CQM) reveals the nature of atoms 
and molecules using classical physical laws for the first time. As discussed above, traditional 
quantum mechanics can solve neither multi-electron atoms nor molecules exactly. By 
contrast, CQM produces exact, closed-form solutions containing physical constants only for 
even the most complex atoms and molecules. 

The present invention is the first and only molecular modeling program ever built on 
the CQM framework. All the major functional groups that make up most organic molecules 
have been solved exactly in closed-form solutions with CQM. By using these functional 
groups as building blocks, or independent units, a potentially infinite number of organic 
molecules can be solved. As a result, the present invention can be used to visualize the exact 
3D structure and calculate the heat of formation of almost any organic molecule. 

For the first time, the significant building-block molecules of chemistry have been 
successfully solved using classical physical laws in exact closed-form equations having 
fundamental constants only. The major functional groups have been solved from which 
molecules of infinite length can be solved almost instantly with a computer program. The 
predictions are accurate within experimental error for over 375 exemplary molecules. 

Applicant's CQM is the theory that physical laws (Maxwell's Equations, Newton's 
Laws, Special and General Relativity) must hold on all scales. The theory is based on an often 
overlooked result of Maxwell's Equations, that an extended distribution of charge may, under 
certain conditions, accelerate without radiating. This "condition of no radiation" is invoked to 
solve the physical structure of subatomic particles, atoms, and molecules. 

In exact closed-form equations with physical constants only, solutions to thousands of 
known experimental values arise that were beyond the reach of previous outdated theories. 
These include the electron spin, g-factor, multi-electron atoms, excited states, polyatomic 
molecules, wave-particle duality and the nature of the photon, the masses and families of 
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fundamental particles, and the relationships between fundamental laws of the universe that 
reveal why the universe is accelerating as it expands. CQM is successful to over 85 orders of 
magnitude, from the level of quarks to the cosmos. Applicant now has over 65 peer-reviewed 
journal articles and also books discussing the CQM and supporting experimental evidence. 

The molecular modeling market was estimated to be a two-billion-dollar per year 
industry in 2002, with hundreds of millions of government and industry dollars invested in 
computer algorithms and supercomputer centers. This makes it the largest effort of 
computational chemistry and physics. 

The present invention's advantages over other models includes: Rendering true 
molecular structures; Providing precisely all characteristics, spatial and temporal charge 
distributions an4 energies of every electron in every bond, and of every bonding atom; 
Facilitating the identification of biologically active sites in drugs; and Facilitating drug 
design. 

An objective of the present invention is to solve the charge (mass) and current-density 
functions of polyatomic molecules, polyatomic molecular ions, diatomic molecules, 
molecular radicals, molecular ions, or any portion of these species from first principles. In an 
embodiment, the solution for the polyatomic molecules, polyatomic molecular ions, diatomic 
molecules, molecular radicals, molecular ions, or any portion of these species is derived from 
Maxwell's equations invoking the constraint that the bound electron before excitation does 
not radiate even though it undergoes acceleratioa 

Another objective of the present invention is to generate a readout, display, or image 
of the solutions so that the nature of polyatomic molecules, polyatomic molecular ions, 
diatomic molecules, molecular radicals, molecular ions, or any portion of these species be 
better understood and potentially applied to predict reactivity and physical and optical 
properties. 

Another objective of the present invention is to apply the methods and systems of 
solving the nature of polyatomic molecules, polyatomic molecular ions, diatomic molecules, 
molecular radicals, molecular ions, or any portion of these species and their rendering to 
numerical or graphical form to all atoms and atomic ions. 

These objectives and other objectives are obtained by a system of computing and 
rendering the nature of at least one specie selected from a group of diatomic molecules 
having at least one atom is other than hydrogen, polyatomic molecules, molecular ions, 
polyatomic molecular ions, or molecular radicals, or any functional group therein, comprising 
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physical, MaxweUian solutions of charge, mass, and current density functions of said specie, 
said system comprising processing means for processing physical, Maxwellian equations 
representing charge, mass, and current density functions of said specie; and an output device 
in communication with the processing means for displaying said physical, Maxwellian 
solutions of charge, mass, and current density functions of said specie. 

Also provided is a composition of matter comprising a plurality of atoms, the 
improvement comprising a novel property or use discovered hy calculation of at least one of 
a bond distance between two of the atoms, a bond angle between three of the atoms, and 
a bond energy between two of the atoms, orbital intercept distances and angles,charge-density 
functions of atomic, hybridized, and molecular orbitals, the bond distance, bond angle, and 
bond energy being calculated from physical solutions of the charge, mass, and current density 
functions of atoms and atomic ions, which solutions are derived from Maxwell's equations 
using a constraint that a bound electrons) does not radiate under acceleration. 

The presented exact physical solutions for known species of the group of polyatomic 
molecules, polyatomic molecular ions, diatomic molecules, molecular radicals, molecular 
ions, or any functional group therein, can be applied to other species. These solutions can be 
used to predict the properties of other species and engineer compositions of matter in a 
manner which is not possible using past quantum mechanical techniques. The molecular 
solutions can be used to design synthetic pathways and predict product yields based on 
equilibrium constants calculated from the heats of formation. Not only can new stable 
compositions of matter be predicted, but now the structures of combinatorial chemistry 
reactions can be predicted. 

Pharmaceutical applications include die ability to graphically or computationally 
render the structures of drugs that permit the identification of the biologically active parts of 
the specie to be identified from the common spatial charge-density functions of a series of 
active species. Novel drugs can now be designed according to geometrical parameters and 
bonding interactions with the data of the structure of the active site of the drug. 

The system can be used to calculate conformations, folding, and physical properties, 
and the exact solutions of the charge distributions in any given specie are used to calculate the 
fields. From the fields, the interactions between groups of the same specie or between groups 
on different species are calculated wherein the interactions are distance and relative 
orientation dependent The fields and interactions can be determined using a finite-element- 
analysis approach of Maxwell's equations. 
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Embodiments of the system for performing computing and rendering of the nature of 
the polyatomic molecules, polyatomic molecular ions, diatomic molecules, molecular 
radicals, molecular ions, or any portion of these species using the physical solutions may 
comprise a general purpose computer. Such a general purpose computer may have any 
number of basic configurations. For example, such a general purpose computer may 
comprise a central processing unit (CPU), one or more specialized processors, system 
memory, a mass storage device such as a magnetic disk, an optical disk, or other storage 
device, an input means, such as a keyboard or mouse, a display device, and a printer or other 
output device. A system implementing the present invention can also comprise a special 
purpose computer or other hardware system and all should be included within its scope. 

Although not preferred, any of the calculated and measured values and constants 
recited in the equations herein can be adjusted, for example, up to + 10%, if desired. 

9 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 illustrates a block diagram of an exemplaiy software program; and 
Figs. 2 and 3 illustrate pictures of an exemplary software program. 
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DETAILED DESCRIPTION OF DSfVETIQN 

Organic molecules comprising an arbitrary number of atoms can be solved using the 
similar principles and procedures as those used to solve alkanes of arbitrary length. Alkanes 
can be considered to be comprised of the functional groups of CH X , CK , and C - C . TTiese 
groups with the corresponding geometrical parameters and energies can be added as a linear 
sum to give the solution of any straight chain alkane as shown in the Continuous-Chain 
Alkanes section. Similarly, the geometrical parameters and energies of all functional groups 
such as alkanes, branched alkanes, alkenes, branched alkenes, alkynes, alkyl fluorides, alkyl 
chlorides, alkyl bromides, alkyl iodides, alkene halides, primary alcohols, secondary alcohols, 
tertiary alcohols, ethers, primary amines, secondary amines, tertiary amines, aldehydes, 
ketones, carboxylic acids, carboxylic esters, amides, N-alkyl amides, N,N-dialkyl amides, 
urea, acid halides, acid, anhydrides, nitriles, thiols, sulfides, disulfides, sulfoxides, sulfones, 

* ' ' - . » ' ' ' 

sulfites, sulfates, nitro aJkanes, nitrites, nitrates, conjugated polyenes, aromatics, heterocyclic 
aromatics, substituted aromatics, and others can bp solved. The fimctional-group solutions 
can be made into a linear superposition and sum, respectively, to give the solution of any 
organic molecule. The solutions of the functional groups can be conveniently obtained by 
using generalized forms of the geometrical and energy equations. The equations and sections 
that are referenced by not contained in this text refer to those sections and equations of the 
book by R. L. Mills entitled, 'The Grand Unified Theory of Classical Quantum Mechanics", 
June 2006 Edition, Cadmus Professional Communications-Science' Press Division, Ephrata, 
PA, ISBN 0963517171, Library of Congress Control Number 2005936834; posted at 

http://w w.blacklightpower.coni/bookdownlQad^html which is incorporated in its 

entirety by reference. ) 

. . - . . ■ . . . 

Consider the case wherein at least two atomic orbital hybridize as a linear 
combination of electrons at the same energy in order to achieve a bond at an energy 
minimum, and the sharing of electrons between two or more such orbitals to form a MO 
permits the participating hybridized orbitals to decrease in energy through a decrease in the 
radius of one or more of the participating orbitals. The force generalized constant k* of a 
H 2 -type ellipsoidal MO due to the equivalent of two point charges of at the foci is given by: 
■ C t (^^ 2c , . . 

'--sir • <«•» 
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where C, is the fraction of the # 2 -type ellipsoidal MO basis function of a chemical bond of 
the molecule or molecular ion which is 0.75 (Eq. (13,59)) in the case of H bonding to a 
central atom and 0.5 (Eq. (14.152)) otherwise, and C 2 is the factor that results in an 

equipotential energy match of the participating at least two molecular or atomic orbitals of the 
chemical bond. From Eqs. (13.58-13.63), the distance from the origin of the MO to each 
focus c 1 is given by: 



c' = a\ ^ -J-iS- (15.2) 
'm/TC&a \2Cf 2 

The internuclear distance is 



2C--2J-23L. (15-3) 
2Cf 2 



The length of the semiminor axis of the prolate spheroidal MO b = c is given by 

b=te-S 2 (15.4) 
And, the eccentricity, e , is 



e = — 
a 



(15.5) 



From Eqs. (1 1.207-1 1.212), the potential energy of the two electrons in the central field of the 
nuclei at the foci is 



~f (15.6) 

O 



%7tS 
c 

The potential energy of the two nuclei is 



2 



?— < 15/7) 

WVa^-6 2 

The kinetic energy of the electrons is 

ft 2 . aWa 2 -£ 2 . 
r = njC,c 2 1 in , (15.8) 

2m as a 2 - 6 2 a - Vtf 2 - 6 2 

* * 

And, the energy j J^, , of the magnetic force between the electrons is 

F » = w*: — n^rr 1 * — 7T^ (15 - 9) 

4roaVa 2 -& 2 a-y]cr~b 2 
The total energy of the if.-type prolate spheroidal MO, £ r is given by the sum of the 
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energy terms: 

E T faMo)=V +T + V +V 

T V * / c m p 



(15.10) 



rifi' 



c^l 2—*- 



\ a+Jcf-b 2 
In- 



a-4a 2 -b 2 



-1 



C ! C 2 



(15.11) 



wbere w, is the number of equivalent bonds of the MO and applies in the case of functional 
groups. In the case of independent MOs not in contact with the bonding atoms, the terms 
based on charge are multiplied by c B0 , the bond-order factor. It is 1 for a single bond, 4 for 

■ » 

an independent double bond as in the case of the C0 2 and N0 2 molecules, and 9 for an 

* 

independent triplet bond. Then, the kinetic energy term is multiplied by c' B0 which is 1 for a 
single bond, 2 for a double bond, and 9/2 for a triple bond c, is the fraction of the /? 2 -type 

ellipsoidal MO basis function of an MO which is 0.75 (Eqs. (13.67-13.73)) in the case of H 

< ■ . 

bonding to an unhybridized central atom and 1 otherwise, and c 2 is the factor that results in 

an equipotential energy match of the participating the MO and the at least two atomic orbitals 
of the chemical bond. Specifically, to meet the equipotential condition and energy matching 

conditions for the union of the H t -type-ellipsqidal-MO and the HOs or AOs of the bonding 

• * > •• . 

atoms, the factor c 2 of a H 2 -rtype ellipsoidal. MO may given by CO . one, (ii) the ratio of the 

Coulombic or valence energy of. the AO or HO of at least one atom of the bond and 
13.605804 eV , the Coulombic energy between the electron and proton of H , (iii) the ratio of 
the valence energy of the AO or HO of one atom and the Coulombic energy of another, (iv) 

. * * * * * * 

i i ■ i • 

the ratio of the valence energies of the AOs or HOs of two atoms, (v) the ratio of two c 2 

factors corresponding to any of cases (ii)-(iv), and (vi) the product of two different c 2 factors 

,'•*■'*••• ■ . ■ 

corresponding to any of the cases (i>(v). Specific examples of the factor c 2 of a # 2 -type 

* * • 

ellipsoidal MO given in previous sections are 

0.936127 , the ratio of the ionization energy of N 14.53414 eV and 13.605804 eV , the 
Coulombic energy between the electron and proton of H ; 

0.91771 , the ratio of 14.82575 eV, -E Coulomb (c,2sp 2 ), and 13.605804 eV ; 
. 0.87495 , the ratio of 15.55033 eV , ~E c ^(£ etJm >2sp*y and 13.605804 eV ; 
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0.85252 , the ratio of 15.95955 eV , -£ CoutoBi (c alvlem ,2sp 3 ), and 1 3.605804 eV ; 
0.85252 , the ratio of 15.95955 eV, -E Caulamt (C bemm ,2sp i \ and 13.605804 eV, and 

0.86359 , the ratio of 15.55033 eV, -E Cmlmt (c ' allam ,2sp 1 ') > and 13.605804 eV . 

In the generalization of the hybridization of at least two atomic-orbital shells to form a 
shell of hybrid orbitals, the hybridized shell comprises a linear combination of the electrons 
of the atomic-orbital shells. The radius of the hybridized shell is calculated from the total 
Coulombic energy equation by considering that the central field decreases by an integer for 
each successive electron of the shell and that the total energy of the shell is equal to the total 

Coulombic energy of the initial AO electrons. The total energy E T (atom 9 msp*^ (m is the 

integer of the valence shell) of the AO electrons and the hybridized shell is given by the sum 
of energies of successive ions of the atom over the n electrons comprising total electrons of 

■ i 

the at least one AO shell. 

■ 

E T (atom,msp*\- (15.12) 

m=l 

where IP m is the wth ionization energy (positive) of the atom. The radius of the 
hybridized shell is given by: 

v - £ r Tc'f >\ 1 • (15 - I3) 

^ q-z- n 8^ 0 E T [atom, msp ) 

• .. _ ■ ■ 

Then, the Coulombic energy E^ omb (atom 9 rnsp 3 J of the outer electron of the atom msp* 
shell is given by 



.2 



E co^(e tom ^ pi y (1514) 

» 

» ■ 

* • * 

In the case that during hybridization at least one of the spin-paired AO electrons is unpaired 
in the hybridized orbital (HQ), the energy change for the promotion to the unpaired state is 
the magnetic energy E(magnetic) at the initial radius r of the AO electron: 

E(magnetic) = ^^J^A (15 .15) 

mr r 

e 

Then, the energy E(atom,msp 3> ) of the outer electron of the atom msp 3 shell is given by the 
sum of E Coulomb (atom, msp 3 ) and E(magnetic): 
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2 2 2 

E(atom,msp 3 \=^—^^fi- (15.16) 

Consider next that the at least two atomic orbitals hybridize as a linear combination of 
electrons at the same energy in order to achieve a bond at an energy minimum with another 
atomic orbital or hybridized orbital. As a further generalization of the basis of the stability of 
the MO, the sharing of electrons between two or more such hybridized oAitals to form a MO 
permits the participating hybridized orbitals to decrease in energy through a decrease in the 
radius of one or more of the participating orbitals. In this case, the total energy of the 

hybridized orbitals is given by the sum of E(atom 9 msp^ and the next energies of successive 

ions of the atom over the n electrons comprising the total electrons of the at least two initial 

AO shells. Here, E(atom,msp 3 ^ is the sum of the first ionization energy of the atom and the 

hybridization energy. An example of E(atom 5 msp 3 ^ for JS(c,2ap 3 ) is given in Eq. (14.503) 
where the sum of the negative, of the first ionization energy of C, -11.27671 eV , plus the 
hybridization energy to form the C2sp* shell given by Eq. (14.146) is 

* 

E(c,2sp 3 )= -14163489 eV . 

>.■.'■••• ■ ■ ■ 

■I t ^ 

Thus, the sharing of electrons between two atom msp 3 HOs to form an atom-atom- 

bond MO permits each participating hybridized orbital to decrease in radius and energy. In 

order to further satisfy the potential, kinetic, and orbital energy relationships, each atom msp 

HO donates an excess of 25% per bond of its electron density to the atom-atom-bond MO to 
form an energy minimum wherein the atom-atom bond comprises one of a single, double, or 

> triple bond. In each case, the radius of the hybridized shell is calculated from the Coulombic 

» 

energy equation by considering that the central field decreases by an integer for each 
successive electron of the shell and the total energy of the shell is equal to the total 
Coulombic energy of the initial AO electrons plus the hybridization energy. The total energy 

E T (moljatom,msp^ {m is the integer of the valence shell) of the HO electrons is given by 

> the sum of energies of successive ions of the atom over the n electrons comprising total 
electrons of the at least one initial AO shell and the hybridization energy: 

E T [moLatom,msp z y^ E{atom 9 msp 3 y £ ZP (15.17) 

m-1 
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where lP m is the mih ionization energy (positive) of the atom and the sum of -IP } plus the 

hybridization energy is E^tom^msp 3 ^. Thus, the radius of the hybridized shell is given 
by: 

r V = 2 (Z-q)s(0.25)\ - , e ^ (15.18) 

\q=z-n ) X7re 0 E r [moLatom, msp \ 

where s = 1,2,3 for a single, double, and triple bond, respectively. The Coulombic energy 

^coulomb (*nol<&tom,msp 3 ^ of the outer electron of the atom msp 3 shell is given by 

E o*m* (molAtom,msp 3 y (15.19) 

0 msp* 

In the case that during hybridization at least one of the spin-paired AO electrons is unpaired 
in the hybridized orbital (HO), the energy change for the promotion to the unpaired state is 
the magnetic energy E(magnetic) at the initial radius rof the AO electron given by Eq. 

(15.15). Then, the energy E^olMtorr^msp 3 ^ of the outer electron of the atom msp 3 shell is 

■ * • 

given by the sum of E^^ (molatom.msp^ mA E{magneiic) : 

E (mohatom, msp 3 )= — — — + 2 ^f * (15 .20) 

E T (atom- atom^msp-y.iht energy change of each atom msp 3 shell with the formation of the 

atom-atom-bond MO is given by the difference between E(molxaom 9 msp 3> ] and 

« 

E(atom,msp 3 y 

E T (atom- atom,msp 3 y= E(molMtom,msp 3 ^ E^tom^msp^ (15.21) 
As examples from prior sections, E CaM (molatom,msp 3 ) is one of: 

■ 

E Coulomb \fathme 9 )> 

. . E c^ faom, ™p' ) is one of E M (c, 2sp> ) and £ CooW (ci, V ); 
, E(mol. a tom,msp 3 ) is .one. of ^(c^,^ 3 ) E^^sp 3 ), 



* < 
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E(atom 9 msp^ is one of and 

E T (atom -atom,msp^ is one of £(c~C,2$p 3 ), £(C = C,2jp 3 ), and. S(Cs C,2sp 3 ); 

atom msp 3 is one of C2sp 3 , C/3^p 3 
£ r (atom - atom (i, ^ insp 5 ) is £ r (C-C,2sp 3 ) and E T (atom-atom(s^msp 3 ^ is 

£ r (c = C,25p 3 ),and 

In the case of the Clsp 3 HO, the initial parameters (Eqs. (14.142-14.146)) are 



^ (Z-n)e l \0e l AA<( „ ( 

; , = X <: = 7 * = 0.91771cr 0 (15.22) 

^8«s 0 (el48.25751eK) ine (el48.25751 eV) 0 



.2 -.2 



£ rmfc . (C,2sp 3 )= — = - = -14.82575 eV (15.23) 

^ / 8?rf 0 r v 8^ 0 0.91771a 0 

.E{magnetic)^ 2 ^° e *f = 8;r ^ = 0.19086 (15.24) 

m 2 (r 3 J (o.84317aJ 



E(c.,2sp 3 y 



-e 2 + 2 W 2 h 2 



= -14.82575 eK + 0.19086 , (15.25) 

= -14.63489 eK 

InEq. (15.18), • 



z-i 



£(Z-?) = 10 (15.26) 

Eqs. (14.147) and (15.17) give 

E T (mol.atom,msp i y= E^^^sp 3 ^ -151.61569 eV (15.27) 

Using Eqs. (15.18-15.28), the final values of , £ Co ^(C2#? 3 ), and E^lsp 2 ^ and the 
resulting E^C-CCTsp^j of the MO due to charge donation from the HO to the MO where 
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BO 



C-C refers to the bond order of the carbon-carbon bond for different values of the 
parameter s are given in Table 15.1. 



Table 15.1. The final values of r Q ^, E 0mM (c2sp^y and J?(C2ap 3 ) and the resulting 

E r ^c - C,C2sp'j of the MO due to charge donation from the HO to the MO where C - C 
refers to the bond order of the carbon-carbon bond. 



MO 
Bond 
Order 
(BO) 


s 1 


s 2 


r cv( a «) 
Final 

> 


(eV) 
Final 


£(C2j5p j ) 

(eV) 
Final 


( BO \ 

E T \C - C,C2sp } j 

(eV) 


I 


1 


0 


0.87495 


-15.55033 


-15.35946 


-0.72457 


n 


2 


0 


0.85252 


-15.95955 


-15.76868 


-1.13379 


m 


3 


0 


0.83008 


. -16.39089 


-16.20002 


-1.56513 


IV 


4 


0 


0.80765 


-16.84619 


-16.65532 


-2.02043 . 



In anpther generalized case of the basis of forming a minimum-energy bond with the 
constraint that it must meet the energy matching condition for all MOs at all HOs or AOs, the 

energy E(molAt6m, msp 3 ^ of the outer electron of the atom msp 3 shell of each bonding atom 

« 

must be the average of E.(moLatom 9 msp 3 ^ for two different values of s : 



E(mol. 



^V^^ J ^^' ^P, ) • (15.28); 



♦ * 

In this case, E T (atom-atom 9 msp 3 ^ the energy change of each atom msp 3 shell with the 

■ 

formation of each atom-atom-bond MO, is average for two different values of s : 

■ 

E T (atom - atomfa ^ msp* E T (atom - atom (s 2 ^ msp 3 ) 



E T (atom - atom, msp 1 ^= -r 



(15.29) 



Consider an aromatic molecule such as benzene givbn in the Benzene Molecule 
section. Each C = C double bond comprises a linear combination of a factor of 0.75 of four 

paired electrons (three electrons) from two sets of two C2sp 3 HQs of the participating carbon 

atoms. Each C-H bond of CH having two spin-paired electrons, one from an initially 



18 



WO 2007/053486 



PCT/US2006/042140 



impaired electron of the carbon atom and the other from the hydrogen atom, comprises the 
linear combination of 75% H 2 -type ellipsoidal MO and 25% C2sp 3 HO as given by Eq. 
(13.439). However, E r (atom- atom,msp^ of the C-#-bond MO is given by 

0.5£ r (c = C,2sp*^ (Eq, (14.247)) corresponding to one half of a double bond that matches 

the condition for a single-bond order for C - H that is lowered in energy due to the aromatic 
character of the bond. 

A further general possibility is that a minimum-energy bond is achieved with 
satisfaction of the potential, kinetic, and orbital energy relationships by the formation of an 
MO comprising an allowed multiple of a linear combination of 7/ 2 -type ellipsoidal MOs and 
corresponding HOs or AOs that contribute a corresponding allowed multiple (e.g. 0.5, 0.75, 
1) of the bond order given in Table 15.1. For example, the alkane MO given in the 

Continuous-Chain Alkanes section comprises a linear combination of factors of 0.5 of a 

■ 

single bond and 0.5 of a double bond. 

; Consider a first MO and its HOs comprising a linear combination of bond orders and 
a second MO that shares a HO with the first. In addition to the mutual HO, the second MO 
comprises another AO or HO having a single bond order or a mixed bond order. Then, in 

order for the two MOs to be energy matched, the bond order of the second MO and its HOs or 

• • • i 

its HO and AO is a linear combination of the terms corresponding to the bond order of the 

* i * . ■ ■ 

mutual HO and the bond order of the independent HO or AO. Then, in general, 

***** • ■ • ■ f ^ • i « 

E T [atom-atom,msp 5 \ the energy change of each atom msp 3 shell with the formation of 

each . atom-atom-bond MO, is a weighted linear sum for different values of s that matches the 
energy of the bonded MOs, HOs, and AOs: 



where c_ is the multiple of the BO of s„ . The radius r of the atom msp 3 shell of each 

■ 

bonding atom is given by the Coulombic energy using the initial energy E Coahnib (atom,msp 3 } 

and E T (atom-atom^msp^ the energy change of each atom msp 3 shell with the formation 
of each atom-atom-bond MO: 

~ e "* 
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where E Cndttmb (C2jp 3 )= -14.825751 eV . The Coulombic energy E CoM (moLatom 9 msp 3 ^ 

of the outer electron of the atom msp 2 shell is given by Eq. (15.19). In the case that during 
hybridization, at least one of the spin-paired AO electrons is unpaired in the hybridized 
orbital (HO), the energy change for the promotion to the unpaired state is the magnetic energy 
E(magnetic) (Eq. (15.15)) at the initial radius rof the AO electron. Then, the energy 

E (moLatom,Msp*^ of the outer electron of the atom msp* shell is given by the sum of 
E coulomb (molMtom, msp* ) and E(magnetic) (Eq. (15.20)). E T (atom — atom,msp 3 ^, the 
energy change of each atom msp 3 shell with the formation of the atom-atom-bond MO is 
given by the difference between E^molMtom.msp^ and E(atom,msp 3 ^ given by Eq. 

► * * ■ 

(15.21). Using Eq. (15.23) for E Cailmb (c,2sp 3 ) in Eq. (15.31), the single bond order energies 

given by Eqs. (15.18-15.27) and shown in Table 15.1, and the linear combination energies 
(Eqs. (15.28-15.30)), the parameters of linear combinations of bond orders and linear 

■ 

combinations of mixed bond orders are given in Table 1 5.2, 

Table 152. The final values of r v , ^^(p*/? 3 ), and E(cisp^ and the resulting 

* j ■ - 

'• BO \ 

& T \C-C,C2sp J of the MO comprising a linear combination of i/ 2 -type ellipsoidal MOs and 

corresponding HOs of single or mixed bond order where c Sn is the multiple of the bond order 
parameter E T (atom - atom (s n \ msp 3 given in Table 15.1. 
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MO 
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^ z 
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1 /IT 
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a ^ 

U.J 


A 
U 
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1 /OTT 

1/Z11 
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A ^ 


a 
u 


n 

u 






1 


U.3 


z 




A 


17211 + (1 + 11) 


z 


A 

U.2D 


1 
1 


A 


z 


^ //ITT 

3/4U 




A *7< 


u 


U 


A 
U 


1 T 11 


i 
1 


U. J 


z 


ft 5 


A 


i+m 


1 


0.5 


3 


0.5 


0 


i+rv 

* 


1 


0.5 


4 


0.5 


0 


n+m 


2 


0.5 


3 


0.5 


0 


n + iy 


2 


0.5 


4 


0.5 


0 


rn+iv 

i 


3 


0.5 


4 


0.5 


0 


IV + IV 


4 


0.5 


4 


0.5 


0 
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MO 
Bond Order 
(BO) 


C 
*> 


C2sp 3 V 0/ 

Final 


E (C2sd 3 ^ 

Coulomb \C ^ J 

(eV) 
Final 


(eV) 
Final 


( BO \ 

. £ r i C-C f C2sp 3 j 
(eV) 


1 y^>f 
1/21 


a 

0 


A OA£*00 

0.89582 


-15.18804 


-14.99717 


-0.36228 


1/211 


A 

0 


0.88392 


-15.39265 


-15.20178 


-0.56689 


1 1 1 /ITT 

1 + 1/211 


0 


0.87941 


-15.47149 


-15.28062 


-0.64573 


1 /OTT _f_ /7 i rrv 
1/2II+ (1 + 11) 


0.25 


0.87363 


-15.57379 


-15.38293 


-0.74804 


^ //ITT 

3/411 


0 


0.86793 


-15.67610 


-15.48523 


-0.85034 


T -U TT 
1 i u 








-I J.D04U/ 


A CsTi A 1 o 

-0.92918 


i+m 


0 


0.85193 


-15.97060 


-15.77974 


-1.14485 


I + IV 


0 


0.83995 


-16.19826 


-16.00739 


-1.37250 


n+m 


0 


0.84115 


-16.17521 


-15.98435 


-1.34946 


n+rv 


0 


0.82948 . 


-16.40286 


-16.21200 


-1.57711 


m+iv 


0 


0.81871 


-16.61853 


-16.42767 


-1.79278 


IV + IV 

• 

> 


0 


0.80765 


-16.84619 
_ 


-16.65532 


-2.02043 



. Consider next the radius of the AO or HO due to the contribution of charge to more 

than one bond. The energy contribution due to the charge donation at each atom such as 

• ' ■ , • ■ ■ 

carbon superimposes linearly,. In general, the radius* of the C2sp 3 HO of a carbon 

atom of a given molecule is calculated using Eq. (14.514) by considering ^E T (MO.lsp 3 ^ 

the total energy donation to each bond with which it participates in bonding, the general 
equation for the radius is. given by 



-e 2 



moi2sp 2 



2 (15.32) 



e 



8^ 0 (el4.82575UK + X|^(MO,2^ 3 )) 

The Coulombic energy ' E Coutomb (molatom 3 msp^ of the outer electron of the atom msp 3 shell 
is given by Eq. (15.19). In the case that during hybridization, at least one of the spin-paired 
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AO electrons is unpaired in the hybridized orbital (HO), the energy change for the promotion 
to the unpaired state is the magnetic energy E(magneti'c) (Eq. (15.15)) at the initial radius 

- • 

rof the AO electroa Then, the energy E(moLatom,msp*) of the outer electron of the 

atommsp 3 shell is given by the sum of E CmlmA (mol.atom,msp i ^ and E(magnetic) (Eq. 
(15.20)). 

For example, the C2sp 3 HO of each methyl group of an alkane contributes 
-0.92918 eV (Eq. (14.513)) to the conesponding single C-C bond; thus, the 
corresponding C2sp 3 HO radius is given by Eq. (14.514). The Clsp 1 HO of each methylene 
group of C a H 2n+2 contributes -0.92918 eV to each of the two corresponding C-C bond 

MOs: Thus, the radius (Eq. (15.32)), the Coulombic energy (Eq. (15.19)), and the energy 
(Eq. (15.20)) of each alkane methylene group are 



-e 2 



— W-y ^o^co^^yZ^T (methylene C-C. V)) 



e 2 



8^e 0 (el4.825751 eV + e0.92918 eV + e0.9291Z eV) 
= 0.81549« 0 



(15.33) 



?<^(c_ w -.V> 5 -^-, = -16.68412 8 K .. (1534) 



• • 2 2''* 

^^ iT 7K^To \ + H ^ =^16.49325 eV (15,35) 

V ' 8as 0 (0.81549fl 0 ) «£(&.84317flJ 

In the determination of the parameters of functional groups, heteroatoms bonding to 
Clsp 3 HOs to form MOs are energy matched to the Clsp* HOs. Thus, the radius and the 
energy parameters of a bonding heteroatom are given by the same equations as those for 
Clsp 2 HOs. Using Eqs. (15.15), (15.19-15.20), (15.24), and (15.32) in a generalized fashion, 

the final values of the radius of the HO or AO, r AtomHOAO , E Cmhmb (moLatom 9 msp 3 y and 
E (pmoi 2s P l J me calculated using £ E T (MO, 2sp 2 \ the total energy donation to each bond 
with which an atom participates in bonding corresponding to the values of £ \C-C s C2sp*\ °f 
the MO due to charge donation from the AO or HO to the MO given in Tables 15.1 and 15.2. 



23 



WO 2007/053486 



PCT/US2006/042140 



Table 15.3.A. The final values of 

r Atom.HOAO 9 ^Coulomb (j W ^ ,fi ^ WW > TrtS P 

calculated using the values of e t \C-C,C2si?j given in Tables 15.1 and 15.2. 
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Table 15.33. The final values of r JbmJguo9 E^^fmoljitom.msp 3 ^ and ^(c^Zsp 3 ) 
calculated for heterocyclic groups using the values of E r ^c*- c,C2sp*j given in Tables 15,1 



and 15.2. 
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The energy of the MO is matched to each of the participating outermost atomic or 
hybridized orbitals of the bonding atoms wherein the energy match includes the energy 
contribution due to the AO or HO's donation of charge to the MO. The force constant k 1 
(Eq. (15.1)) is used to determine the ellipsoidal parameter c % (Eq. (15.2)) of the each H 7 - 
type-ellipsoidal-MO in terms of the central force of the foci. Then, c x is substituted into the 
energy equation (from Eq. (15.11))) which is set equal to «, times the total energy of H 2 

where is the number of equivalent bonds of the MO and the energy of H 2 , 

-31.63536831 eV 9 Eq. (11.212) is the minimum energy possible for a prolate spheroidal 
MO. From the energy equation and the relationship between the axes, the dimensions of the 
MO are solved. The energy equation has the semimajor axis a as it only parameter. The 
solution of the semimajor axis a then allows for the solution of the other axes of each prolate 
spheroid and eccentricity of each MO (Eqs. (15.3-15.5)). The parameter solutions then allow 
for the component and total energies of the MO to be determined. 

The total energy, E T (p 2 Mo), is given by the sum of the energy terms (Eqs. (15.6- 
15.11)) plus E T (AO/HO): 

E T (p 2 Mo)=V e + T + V m + V p + E r (AO/HO) ' (15.36) 



E T (h 2 mo)= 



C l C 7 



e t (ao/ho) 



C l C 2 



2-^ In 



a) a-c' 



-1 



(15.37) 



+ 



e t (ao/ho) 



where «, is the number of equivalent bonds of the MO, c, is the fraction of the H 2 -type 
ellipsoidal MO basis function of a chemical bond of the group, c 2 is the factor that results in 
an equipotential energy match of the participating at least two atomic orbitals of each 
chemical bond, and E f (AO I H<3) is the total energy comprising the difference of the energy 

E(AO/ HQ) of at least one atomic or hybrid orbital to which the MO is energy matched and 

any enerjgy component AE ff M0 (AO/ HO) due to the AO or HO's charge donation to the 



MO. 



E T (AO I HOy E(AO I HOy AE„ iM0 (AOJ HO) 



(15.38) 
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As specific examples given in previous sections, E r (AO I HO) is one from the group of 
E T (AO / HO)= E(02p shell)= -E(ionization; O) = -13.6181 eV ; 
E T (AOIHO)= E(N2p shell^ -E(ionizatiorr, N) = -14.53414 eV ; 
E T (AO/HO)= E(c,2sp 3 y= -14.63489 eV ; 
E T (AO I HO)= E Coutomb (C/, V)= -14.60295 eV ; 
E T (AO/HO)= E{iomzatiorr, C) + E(ionkation; C + ); 
E T (AO I HO)= E^^y -15.35946 eV ; 
EriAOlHOy+E^^lspyE^,^), 



E T (AO I HO)= E(c,2sp 3 )-2E T (c = C,2ap 3 )= -14.63489 eV - (-2.26758 eV); 

i 

E T (AO I HOy E(£ mvkm , V> E(C a ^,2^ i?(c^ w V> 1620002 eV ; 
E T (AO I HOy E(c,2sp 3 y2E T (c = C,2*/> 3 )= -14.63489 eV - (-3.13026 eV\ 
E T (AO I HO)= £(C W(K , V> 4 W V) 

£ r (40 / E(c,2sp 3 y E T (C = C, V)= -14.63489 eV - (-1.13379 eK) and 
E T (AO I HOy E(C dlme ,2sp 3 y~ 15.56407 eV . 



To solve the bond parameters and energies, c' = a t 
(15.2)) is substituted into E to give 



ft 2 4ffi£, 



m,e*2Cf 2 a V 2C i C j 



(Eq. 



%7ts Va 2 - 



0 
.2 



-1 



£ r (^O / jto) 



8jze 0 <?' 



f a. 
c.cJ 2--S- 



'V2 



J 



ln^-1 
a-c' 



e t (aoiho) 



(15.39) 



8^ 04 /-^- 
04, 2C,C 2 



c,c 2 1 2 



_2 



In 



"1 






-J 




2C,C 2 



-1 



E T (AO I HO) 



The total energy is set equal to E(basis energies) which in the most general case is given by 
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the sum of a first integer Wj times the total energy of H 2 minus a second integer n 2 times the 
total energy of H , minus a third integer w 3 times the valence energy of e(AO) (e.g. 

E(Ny= -14.53414 eV) where the first integer can he 1,2,3..., and each of the second and 

third integers can be 0,1,2,3.... 

E(basis energies) - n x (-31.63536831 eV)- n 2 (-13.605804 eV)- ^e(a6) (15.40) 

In the case that the MO bonds two atoms other than hydrogen, Efbasis energies) is n, times 

the total energy of H 2 where is the number of equivalent bonds of the MO and the energy 

of H 2 , -31.63536831 eV , Eq. (11.212) is the minimum energy possible for a prolate 
spheroidal MO: 

E(basis energies) = ^ (-31.63536831 eV) (15.41) 

E T (p 2 Mo\ is set equal to E(basis energies), and the semimajor axis a is solved. 
Thus, the semimajor axis a is solved from the equation of the form: 



04 '2CC 



I „ a. 
e,c 2 2— 2- 



In 







\2Cf 2 






[2C,C 2 



-1 



+ £ r (AO I HCi)= E(basis energies) (15.42) 



The distance from the origin of the 7f 2 -type-ellipsoidal-MO to each focus c\ the internuclear 
distance 2c ! , and the length of the semiminor axis of the prolate spheroidal' # 2 -type MO 

2> = c are solved from the semimajor axis a using Eqs. (15.2-15.4). Then, the component 
energies are given by Eqs. (15.6-15.9) and (15.39). 

The total energy of the MO of the functional group, E T (wo), is the sum of the total 

energy of the components comprising the energy contribution of the MO formed between the 

participating atoms and E T (atom - atom, msp 3 JtC>y the change in the energy of the AOs or 

HOs upon forming the bond. From Eqs. (15.39-15.40), E T (po) is 

E T (wo)= E(basis energies) + E T (atom - atom 9 msp*.AO) (15.43) 

During bond formation, the electrons undergo a reentrant oscillatory orbit with 
vibration of the nuclei, and the corresponding energy E^ is the sum of the Doppler, E D , and 

average vibrational kinetic energies, E Kvib : 
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E 



\ 



mc 1 2 VpJ 



(15.44) 



where «, is the number of equivalent bonds of the MO , k is the spring constant of the 
equivalent harmonic oscillator, and fi is the reduced mass. The angular frequency of the 

reentrant oscillation in the transition state corresponding to E D is determined by the force 

between the central field and the electrons in the transition state. The force and its derivative 
are given by 



ao 4xe 0 R 3 



(15.45) 



and 



C, C, e 2 
80 4xs„R 3 



(15.46) 



such that the angular frequency of the oscillation in the transition state is given by 



co = \\- 



1 



m 




'BO 



m 



(15.47) 



where R is the semimajor axis a ox the semiminor axis b depending on the eccentricity of 

the bond that is most representative of the oscillation in the transition state, c B0 is the bond- 

• ...» . ' 

■ 

order factor which is 1 for a single bond and when the MO comprises n x equivalent single 
bonds as in the case of functional groups. c B0 is 4 for an independent double bond as in the 
case of the C0 1 and Jid molecules and 9 for an independent triplet bond C lo is the 

• ■ 

fraction of the H 2 -type ellipsoidal MO basis function of the oscillatory transition state of a 

chemical bond of the group, and C 2o is the factor that results in an equipotential energy match 
of the participating at' least two atomic orbitals of the transition state of the chemical bond. 
Typically, C lo = C, and C 2o = C 2 . The kinetic energy, E K , corresponding to E D is given by 
Planck's equation for functional groups: 



E = fto = ^ 



K 



4jts 0 R 3 
m 



(15.48) 



The Doppler energy of die electrons of the reentrant orbit is 
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2h^ 


4th> 0 /J 3 





m c 

e 



E ose given by the sum of E D and E Kvih is 

/ 



E osc n x (E D + E K¥i ^y n, 



E 



fry 













m c 2 

0 



vlb 



(15.49) 



(15.50) 



2?. of a group having ra, bonds is given by £j.(mo)/ «, such that 

/ r-=— r-\ 

IOC i I iu 



. , 2E K 1 It 

£ r (HO)/ «, 4 / f + - ft J — 



Mc 2 2 \n) 



£ T+ fcm V ) is given by the sum of E T (m>) (Eq. (15.42)) and E m (Eq. (15.5 1)): 



7*+qsc 



OJCC 



n,e 



2 



°M2C I C 2 



1\ 



C 1 C 2 



12 — -Jin \ 1 -1 



+£ r (AO t £ r (atom - atom, msp 3 .AO) 



J 



1 + 



2t* 



Ans R l 
\ 



• me 

8 



l 2 \m 



J 



(15.51) 



= (E(basis energies) + E T (atom - atom,msp 3 .AO 









1 4xsR 3 


2^ 




mc 2 

6 



1 - ft 

'2 Hfi 
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(15.52) 

Hie total energy of the functional group E T (pou P ) is the sum of the total energy of the 

components comprising the energy contribution of the MO formed between the participating 
atoms, E(basis energies) y the change in the energy of the AOs or HOs upon forming the 

bond (E T (atom- atom 9 mp 3 .AO)) 9 the energy of oscillation in the transition state, and the 
change in magnetic energy with bond formation, E . From Eq. (15.52), the total energy of 



the group E T (group) is 



E T $Zrovp)= 



(E(basis energies) + E T (atom - atow,m^.^0^1 + 1 



2^ 


C, C, c 2 


0 




m c 1 

e 



(15.53) 



The change in magnetic energy E mag which arises due to the formation of unpaired electrons 
in the corresponding fragments relative to the bonded group is given by 



2 
B 



(15.54) 



where r 3 is the radius of the atom that reacts to form the bond and c 3 is the number of 



electron pairs. 



E„ (group)— 



• * 



(E(basis energies) + E T (atom - atom 5 msp 3 ^o)j 1 + | 



2h 



i 



4tueR 



m 



mc 



V 



(15.55) 



The total bond energy of the group E D (group) is the negative difference of the total energy of 
the group (Eq. (15.55)) and the total energy of the starting species given by the sum of 
c 4 E mlai (c 4 AO I HO) and* c 5 E Mr!a! (c 5 AO I HO): 
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E ^nmp)= - 



(E(basis energies) +• E T (atom - atow7,w^ 3 ^40^1 + j 



2^ 




| 4ite g R 3 


■ m . 


w c 2 



(15.56) 



In the case of organic molecules, the atoms of the functional groups are energy matched to the 
C2sp % HO such that 

E(AO/HO)= -14.63489 eV 



(15.57) 



For examples of E from previous sections: 

r (0.91771aJ 



14803 eV 



(15.58) 



map 



3 r 3 



p- = c,0.11441eK 

a: 



(15.59) 



2 



3 r 3 



" c s , SaM ° MB - = c 3 0.14185 eK 



(15.60) 



(0.93084aJ 

i . " ' 

la the general case of the solution of an organic functional group, the geometric bond 
parameters are solved from the semimajor axis and the relationships between the parameters 
by first using Eq. (15.42) to arrive at a. . Then, the remaining parameters are determined 
using Eqs. (15.1-15.5). Next, the energies are given by Eqs. (15.52-15.59). To meet the 
equipoteiitial condition for the union of the 7/,-type-ellipsoidal-MO and the HO or AO of the 
atom of a functional group, the factor c 2 of a H 2 -typo ellipsoidal MO in principal Eqs. 
(15.42) and (15.52) may given by 



(i) one: 
c 2 =l 



(15.61) 



• (ii) the ratio that is less than one of 13.605804 eV , the magnitude of the Coulombic 
energy between the electron and proton of H given by Eq. (1.243), and the magnitude of the 
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Coulombic energy of the participating AO or HO of the atom, E CM (MOxttom, »isp 3 ) given 
by Eqs. (15.19) and (15.3 1-15.32). For \E (MaiA (MOjatom,msp i )> 13. 605804 eV: 



.2 



ins^cu 13.605804 eV 

c 2 = JT 3 = r- t ^ (1 5.62) 

\ E c<*hmb [MOuttom,msp 3 ) 

in£ »r A . BAerBsp , 



For | £ cw-i (Maa/ow, W ip 3 }<13.605804 eV: 



c 2 = 



**W A - BM J \E Cailanb (MO.atom,msp 3> \ 

e 2 13.605804 eV (15.63) 



Sne 0 a 0 



(in) the ratio that is less than one of 13.605804 eV , the magnitude of the Coulombic 
energy between the electron and proton of H given by Eq. (1.243), and the magnitude of the 
valence energy, E(yalence), of the participating AO or HO of the atom where E(yalence) is 

the ionization energy or E(MOMtom,msp 3 ^ given by Eqs. (15.20) and (15.31-15.32). For 

\E(yaJence)\ >13.605804 eV : 



e 2 



8>re 0 o 0 13.605804 eV 

c i = i* — = ■ ... (15.64) 



% tie jr. _ « 

0 A-B AorBsp' 



\E(vaIence)\ 



For \E(yalence)\ < 13.605804 eV : 



e 1 



= ^o^BAorBsp 3 _ \E(yalence)\ 

13.605804 eV (15 ' 65) 



g 2 



(iv) the ratio that is less than one of the magnitude of the Coulombic energy of the 
participating AO or HO of a first atom, E Cutgmb (MO.atom,msp 3 ^ given by Eqs. (15.19) and 
(15.31-15.32), and the magnitude of the valence energy, E(yalence), of the participating AO 
or HO of a second atom to which the first is energy matched where E(valertce) is the 
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ionization energy or E(MOjatom,msp z ^ given by Eqs. (15.20) and (15.31-15.32). For 
E couiomb (MOxftom,msp 3 ^> E{valence) : 



C 2 = 



\E(vaIence)\ 



(15.66) 



For 



c 2 = 



^coulomb (MOatom,msp*} 
E ccuhmb (MO. atom, msp 3 j < E(yalence) : 

\E(valence)\ 



(v) the ratio that is less than one of the magnitude of the valence-level energies, 
E n (valence) 9 of the AO or HO of the nth participating atom of two that are energy matched 

where Eiyalence) is the ionization energy or E(MO.atom 9 msp 3 ^ given by Eqs. (15.20) and 

(15.31-15.32): 

E^valence) 



c 2 = 



(15.68) 
E 2 (yalence) 



(vi) the factor that is the ratio of the hybridization factor c 2 (l) of the valence AO or 

HO of a first atom and the hybridization factor c 2 \2 j of the valence AO or HO of a second 

atom to which the first is energy matched where c 2 (w) is given by Eqs. (15.62-15.68); 

alternatively c 2 is the hybridization factor c 2 (l) of the valence AOs or HOs a first pair of 

atoms and the hybridization factor c 2 (2) of the valence AO or HO a third atom or second 
pair to which the first two are energy matched: 



c 2 = 



c,(2) < 15 - 69) 



(vii) the factor that is the product of the hybridization factor c 2 (l) of the valence AO 
or HO of a first atom and the hybridization factor c 2 (2) of the valence AO or HO of a second 
atom to which the first is energy matched where c 2 (n) is given by Eqs. (15.62-15.69); 
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alternatively c % is the hybridization factor c 2 (l) of the valence AOs or HOs a first pair of 

atoms and the hybridization factor c 2 (2) of the valence AO or HO a third atom or second 
pair to which the first two are energy matched: 

c 2 = c 2 (l)c 2 (2) (15.70) 

■ 

The hybridization factor c 2 corresponds to the force constant k (Eqs. (11.65) and (13.58)). 

In the case that the valence or Coulombic energy of the AO or HO is less than 13.605804 eV , 
the magnitude of the Coulombic energy between the electron and proton of H given by Eq. 
(1.243), then C 2 corresponding to k } (Eq. (15.1)) is given by Eqs. (15.62-15.70). 

■ 

Specific examples of the fectors c 2 and C 2 of a # 2 -type ellipsoidal MO of Eq. 
(15.51) given in following sections are 

c^HO ,0 F), H£*p^ m y :%gg> W7W). 0.77087 ; 
C 1 (C2^^»a)= i ^^(C2^ O )=5!|^(0. 91 77.)=0.8 1 3l7 i 

/„ ( fW v -14.53414 eV „„„ ' 

(^{Hto l°N)= > > 7 A = = 0.94627; 

7 £(C,2.y/? 3 ) -15.35946 eF 

« fa ^ *M - g W - ~ 14 53414 gF - 0 93383 • 
A ^ *(C.2*>) -15.56407 ~ °' 93383 ' 

C 2 * H) = ^ = - 10 - 36001 ^ = 0.76144 ; 
2V ' E(H) -13.60580 eV 

C, (C2s/m ,o S y rjS^^m,). 0.64955 , 
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E(0) 



(p to S3sp 3 to C2sp 3 HO)=-^c 2 (C2sp 3 HO) 



-10.36001 eV y ' 
= 1.20632 



c , , s ^y -''f 099 - 0.85045 , 

* v ' E(H) -13.60580 eV 

C, (C2sp 3 HO to S3sp 3 )= E 9 3SP l c 2 (S3sp 3 )= ~ U - 52126 g *V o,85045)= 0.66951 ; 
* v " y ' £(C,2sp 3 ) 2V ' -14.63489 eV x ' 

C 2 (S3sp 3 to O to C2sp 3 HO)= ^j^T <a (f2sp 3 HO) 

-11.52126 eV , x 

(0.91771) ; 



-13.61806 eV 
= 0.77641 



C 2 



(O to N2p to C2sp 3 HO)= J©** (C2sp 3 HO) 



= -13.61806 g ( } 
-14.53414 eF v y ' 
= 0.85987 



™ < Ci(C2sp 3 HO to N) 0.9 140 

e 2 (#2/> to 02/?)= -4- — * : — <= — — — = 1.06727 ; 

c 2 (C2sp 3 HO to O) 0.85395 

/ , \ / , ^ 13.605804 e7 

C 2 (benzeneC2sp 3 HO)= c 2 (benzeneC2sp 3 HO)= = 0.85252 ; 

15.95955 0V 



°2 



(arylC2sp 3 HO to 0)= ?' 0 ' 3 < c 2 {arylClsp 3 HO) 



m -13.61806 «K ( } 
-14.63489 eV y ' 
= 0.79329 



,rt E (N) -14.53414^7 

(H to anlinc N)= / \\ = = 0.92171 ; 

2K } E[C y 2sp 3 ) -15.76868 eV 



Cj (arylC2sp 3 HO to N)- f^ N \\ c i (arylC2sp 3 HO) 

E i C 9 "2>sp ) 



-14.53414 eF /noc „ r . 

= (0.85252) , and 

-14.63489 eV y ' 

= 0.84665 
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C 2 (S3 P to aryl-type C2^HO)- ^4 = =^gi ^ = 0.65700 . 
v y £(C,2.sp 3 ) -15.76868 eF 



MO INTERCEPT ANGLES AND DISTANCES 

Consider the general case of Eqs. (13.84-13.95) wherein the nucleus of a B atom and the 
nucleus of a A atom comprise the foci of each i/ 2 -type ellipsoidal MO of an A -B bond. 
The parameters of the point of intersection of each H 2 -type ellipsoidal MO and the -4 -atom 
AO are determined from the polar equation of the ellipse: 

1 + e 



r = r A 



1 + 6COS0' 



(15.71) 



The radius of the A shell is r A , and the polar radial coordinate of the ellipse and the radius of 
the A shell are equal at the point of intersection such that 



1 + — 



(15.72) 



1 + — COS0' 

a 



The polar angle 9 X at the intersection point is given by 



( 



6* = cos -1 



a 



t C 
1 + — 



(15.73) 



Then, the angle 0, ^ the radial vector of the A AO makes with the internuclear axis is 

0AAo = m °- 0t (15.74) 
The distance from the point of intersection of the orbitals to the internuclear axis must be the 
same for both component* orbitals such that the angle at = 0 HMO between the internuclear 

axis and the point of intersection of each # 2 -type ellipsoidal MO with the A radial vector 

* 

obeys the following relationship: 

r A^ d AAo = b&ia9 H lUO (15.75) 
such that 



e 



H 2 MO 



sin — — 



(15.76) 
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The distance d H m along the intemuclear axis from the origin of # 2 -type ellipsoidal MO to 



the point of intersection of the orbitals is given by 

d* 2 Mo = °™ 0 H lM o (15 - 77) 
The distance d A AO along the interauclear axis from the origin of the A atom to the point of 
intersection of the orbitals is given by 

d =c'-d ( 15 - 78 ) 

°AAO C H 2 MO V 



BOND ANGLES 

Further consider an ACB MO comprising a linear combination of C - A -bond and C - B 
bond MOs where C is the general central atom. A bond is also possible between the A and 
B atoms of the C-A and C-B bonds. Such A-B bonding would decrease the C-A 
and C-B bond strengths since electron density would be shifted from the latter bonds to the 
former bond. Thus, the ZACB bond angle is determined by the condition that the total 
energy of the # 2 -type ellipsoidal MO between the terminal A and B atoms is zero. The 
force constant k' of a # 2 -type ellipsoidal MO due to the equivalent of two point charges of 
at the foci is given by: 

ft»_ C i C 2 2g2 (15.79) 

o 

where C x is the fraction of the H 2 -type ellipsoidal MO baisis function of a chemical bond of 
the molecule which is 0.75 (Eq. (13.59)) for a terminal A - H ( A is H or other atom) and 1 

■ 

otherwise and C, is the factor that results in an equipotential energy match of the 
participating at least two atomic orbitals of the chemical bond and is equal to the 
corresponding factor of Eqs. (15.42) and (15.52). The distance from the origin of the MO to 
each focus c' of the A- B ellipsoidal MO is given by: 



C'—Q. 



\^o_ (15.80) 



m/2C x C 2 a \2Cf. 
The intemuclear distance is 



2c' = 2. 



£3>_. (15.81) 



2C,C 2 
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The length of the semiminor axis of the prolate spheroidal A - B MO b = c is given by Eq. 
(15.4). 

The component energies and the total energy, E T fayo), of the A-B bond are given 

■ 

by the energy equations (Eqs. (11.207-11.212), (11.213-11.217), and (11.239)) of H 2 except 
that the terms based on charge are multiplied by , the bond-order factor which is 1 for a 
single bond and when the MO comprises equivalent single bonds as in the case of 
functional groups. c B0 is 4 for an independent double bond as in the case of the C0 2 and 
N0 2 molecules. The kinetic energy term is multiplied by c\ 0 which is 1 for a single bond, 2 
for a double bond, and 9/2 for a triple bond. The electron energy terms are multiplied by c x , 
the fraction of the H 2 -type ellipsoidal MO basis function of a terminal chemical bond which 
is 0.75 (Eq. (13.233)) for a terminal A-H {A is H or other atom) and 1 otherwise. The 
electron energy terms are further multiplied by c' 2 , the hybridization or energy-matching 

r 

factor that results in an equipotential energy match of the participating at least two atomic 
orbitals of each terminal bond. Furthermore, when A-B comprises atoms other than H , 

E T (atom - atom,msp 3 .AO>y the energy component due to the AO or HO's charge donation to 
the terminal MO, is added to the other energy terms to give E t (p 2 Mo): 



-e 
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cjfA 2c_-c' -2-lln 



1~2 



a + c % 



BO 



BO 



Hln 

) a-c\ 



-1 



,+ E T (atom - atom,msp 3 \AO^ 



(15.82) 

The radiation reaction force in the case of the vibration of A - B in the transition state 
corresponds to the Doppler energy, E D , given by Eq. (11.181) that is dependent on the 
motion of the electrons and the nuclei. The total energy that includes the radiation reaction of 
the A-B, MO is given by the sum of E T (p 2 Mo) (Eq. (15.82)) and. E^- given Eqs. (11.213- 

11.220), (11.231-11.236), and (11.239-11.240). Thus, the total energy E t (a-b) of the 
A-B MO. including the Doppler term is 
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e t (a-b)= 
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(15.83) 

where C Xo is the fraction of the i/ 2 -type ellipsoidal MO basis function of the oscillatory 
transition state of the A- B bond which is 0.75 (Eq. (13.233)) in the case of H bonding to a 
central atom and 1 otherwise, C 2o is the factor that results in an equipotential energy match of 
the participating at least two atomic orbitals of the transition state of the chemical bond, and 



is the reduced mass of the nuclei given by Eq. ( 1 1 . 1 54). To match the boundary 



condition that the total energy of the A - B ellipsoidal MO is zero, E T (a - B) given by Eq. 



(15.83) is set equal to zero. Substitution of Eq. (15.81) into Eq. (15.83) gives 
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(15.84) 

The vibrational energy-term of Eq. (15.84) is determined by the forces between the central 
field and the electrons and those between the nuclei (Eqs. (11.141-11.145)). The electron- 
central-field force and its derivative are given by 



f(?> 



-c 



60 \ne,a 3 



(15.85) 



and 
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The nuclear repulsion force and its derivative are given by 



(15.86) 



f(a + c')= j ;r 



(15.87) 



and 



/'(*+*■)=- 



e 



such that the angular frequency of the oscillation is given by 



(15.88) 
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(15.89) 



Since both terms of E m = E D + E Ktlb are small due to the large values of a and c' , to very 
good approximation, a convenient form of Eq. (15.84) which is evaluated to determine the 
bond angles of functional groups is given by 
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(15.90) 



From the energy relationship given by Eq. (15.90) and the relationship between the axes 
given by Eqs. (15.2-15.5), the dimensions of the A-B MO can be solved. The most 
convenient way to solve Eq. (1 5.90) is by the reiterative technique using a computer. 

A factor c 2 of a given atom in the determination of c $ 2 for calculating the zero of the 
total A-B bond energy is typically given by Eqs. (15.62-15.65). In the case of a H-H 
terminal bond of an alkyl or alkenyl group, c' 2 is typically the ratio of c 2 of Eq. (15.62) for 
the H- H bond which is one and c 2 of the carbon of the corresponding C-H bond: 
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c , = 1 _ 13.605804 eV 

C * c 2 (C2V) ^(c-flcV) ' 
In the case of the detennination of the bond angle of the ACH MO comprising a linear 
combination of C - A -bond and C-H -bond MOs where A and C are general, C is the 
central atom, and c 2 for an atom is given by Eqs. (15.62-15.70), c' 2 of the A- H terminal 

bond is typically the ratio of c z of the A atom for the A - H terminal bond and c 2 of the C 
atom of the corresponding C-H bond: 

c 2 (A(A-H)msp>) 

2 c 2 (C(C-HXmsp 3 ) (15 ^ 

In the case of the determination of the bond angle of the COH MO of an alcohol comprising 
a linear combination of C-0-bond and O- //-bond MOs where C, O, and if are carbon, 

oxygen, and hydrogen, respectively, c 2 of the C-H terminal bond is typically 0.91771 
since the oxygen and hydrogen atoms are at the Coulomb potential of a proton and an electron 
(Eqs. (1.236) and (10.162), respectively) that is energy matched to the Clsp 3 HO. 

In the determination of the hybridization factor d % of Eq. (15.90) from Eqs. (15.62- 

15.70), the Conlombic energy, E Coulomb (MO.atom,msp*y or the energy, E^MOxitom^msp^ 

the radius r AB AorB ^ of the A or B AO or HO of the heteroatom of the A-B terminal bond 

MO such as the Clsp 3 HO of a terminal C-C bond is calculated using Eq. (15.32) by 
considering V E T (MO,2sp 3 \ the total energy donation to each bond with which it 
participates in bonding as it forms the terminal bond. The Cpulombic energy 
E Couhmb (MO.atom 3 msp 3 ^ of the outer electron of the atom msp 3 shell is given by Eq. (15.19). 

In the case that during hybridization, at least one of the spin-paired AO electrons is unpaired 
in the hybridized orbital (HO), the energy change for die promotion to the unpaired state is 
the magnetic energy E(magnetic) (Eq. (15.15)) at the initial radius r of the AO electron. 

Then, the energy E^MOMtom.msp 3 ^ of the outer electron of the atom msp 3 shell is given by 
the sum of E^^^MOxitom.msp^ and E(magnetic) (Eq. (15.20)). 

In the specific case of the terminal bonding of two carbon atoms, the c 2 factor of each 
carbon given by Eq. (15.62) is determined using the Coulombic energy 
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E (C-CC2sp l \ of the outer electron of the C2sp* shell given by Eq. (15.19) with the 

Conlomb \ 7 

radius r , of each Clsp 1 HO of the terminal C-C bond calculated using Eq. (15.32) 

C-C Clxjr 

by considering Yj E t (^O.lsp 1 ^ the total energy donation to each bond with which it 

participates in bonding as it forms the terminal bond including the contribution of the 
methylene energy, 0.92918 eV (Eq. (14.513)), corresponding to the terminal C-C bond. 

* 

The corresponding E T (atotn-atom,msp\AOya Eq. (15.90) is 

E T (C-C C2s/> 3 )= -1.85836 eV . 

In the case that the terminal atoms are carbon or other heteroatoms, the terminal bond 
comprises a linear combination of the HOs or AOs; thus, c' % is the average of the 
hybridization factors of the participating atoms corresponding to the normalized linear sum: 

C 2 = ~( C 2 (? t0m 0+ C 2 2 )) (15 ? 3) 



In the exemplary cases of C - C , O-O, md N-N where C is carbon: 



< 0 2 ^ 
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,. , if , , 13.605804 eF ( 13.605804 eF J j 
. .. " 21^ (A- MAO I Hof E CaJomb (A - Aj^AO I HO)} 

In the exemplary cases of C-N, C-O, and C-S, 



(15.94) 



2 2 



13.605804 



\ Coulomb 



(15.95) 



where C is carbon and c 2 (C fo 5) is the hybridization factor of Eqs. (15.52) and (15.84) that 

■ 

matches the energy of the atom B to that of the atom C in the group. For these cases, the 
corresponding: E (atom- atom,msp 3 .Ao) term in Eq. (15.90) depends on the hybridization 

and bond order of the terminal atoms in the molecule, but typical values matching those used 
in the determination of the bond energies (Eq. (15.56)) are 

E T (c - O C2sp 3 .02 P y -1.44915 eV ; E T (c-0 C2sp\02p)= -1 .65376 eV ; 
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E T (C - N C2sp 3 J/2p)= -1 .44915 eV ; E T (c - S C2sp 3 X2 />)= -0.72457 eV ; 

E T (p-0 02p.02p)= -1.44915 eV; E T (pr ° 02p.02p)= -1.65376 eV; 

E T (N-NN2 /?JV2p)= -1 .4491 5 eV ; E T (N-ON2 p.02 p)= -1 .4491 5 eV ; 

E T (F- F F2p.F2p)=-\M9\5 eV ; E T (ci-Cl Cl3p.Cl3p)= -0.92918 eV; 

E T (Br - Br Br4p.Br4 p)= -0.92918 eV ; E T (l - 1 15p.I5p)= -0.36229 eV \ 

E T (C-FC2sp 3 .F2p)=-Lm36 eV; E T (C-C1 C2sp\Cl3p)= -0.9291 8 eV; 

E T (C - Br C2sp\Br 4 p)= -0.72457 eV; E T (C-I C2sp 3 .I5p\= -0.36228 eV , and 
E T (p-Cl 02p.Cl3p)= -0.9291 8 eV . 

In the case that the terminal bond is X - X where X is a halogen atom, c x is one, 

* * 

and c\ is the average (Eq. (15.93)) of the hybridization factors of the participating halogen 

i 

atoms given by Eqs. (15.62-15.63) where E^ ulmb ^MO.atom,msp 3 ^ is determined using Eq. 

(15.32) and E Cmhmb (MO.atom,msp 3 y= 13.605804 eV for X=I. The factor C, of Eq. 

(15.90) is one for all halogen atoms. The factor C 2 of fluorine is one since it is the only 

halogen wherein the ionization energy is greater than that 13.605804 eV , the magnitude of 
the Coulombic energy between the electron and proton of H given by Eq. (1.243). For each 

of the other halogens, CI 9 Br , and J, C 2 is the hybridization factor of Eq. (15,52) given by 
Eq. (15.70) with c 2 (l) being that of the halogen given by Eq. (15.68) that matches the 
valence energy of X ( (valence)) to that of the C2sp 3 HO (E 2 (valence)^ -14.63489 eV, 
Eq. (15.25)) and to the hybridization of C2sp 3 HO (c 2 (2)= 0.91771, Eq. (13.430)). 

E T {atom-atom,msp 3 .A6) of Eq. (15.90) is the maximum for the participating atoms which 
is -1.44915 eV, -0.92918eF, -0.92918 eV , and -0.33582 eV for F, CI, Br,md I, 

■ « 

respectively. 

Consider the case that the terminal bond is C - X where C is a carbon atom and X 
is a halogen atom. The factors c x and C x of Eq. (15.90) are one for all halogen atoms. For 

X-F 9 c' 2 is the average (Eq. (1 5.95)) of the hybridization factors of the participating carbon 

# 

and F atoms where c 2 for carbon is given by Eq. (15.62) and c 2 for fluorine matched to 
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carbon is given by Eq. (15.70) with c 2 (l) for the fluorine atom given by Eq. (15.68) that 

matches the valence energy of F (E^alence^ -17.42282 eV) to that of the C2sp* HO 

(E 2 (valence)=-U.634&9eV 9 Eq. (15.25)) and to the hybridization of C2sp 3 HO 
(c 2 (2)= 0.91771, Eq. (13.430)). The factor C 2 of fluorine is one since it is the only halogen 
wherein the ionization energy is greater than that 13.605804 eV , the magnitude of the 
Coulombic energy between the electron and proton of H given by Eq. (1 .243). For each of 
the other halogens, Cl 9 Br, and I, c\ is the hybridization factor of the participating carbon 
atom since the halogen atom is energy matched to the carbon atom. C 2 of the terminal-atom 
bond matches that used to determine the energies of the corresponding C - X -bond MO. 
Then, C 2 is the hybridization factor of Eq. (15.52) given by Eq. (15.70) with c 2 (l) for the 

halogen atom given by Eq. (15.68) that matches the valence energy of X (E^alence)) to 

that of the C2sp* HO (E 2 (yalence)= -14.63489 eV , Eq. (15.25)) and to the hybridization of 

Clsp* HO (c 2 (2)= 0.91771, Eq. (13.430)). E r (atom -atom, msp 3 .AO) of Eq. (15.90) is 

♦ » 

the maximum for the participating atoms which is -1.85836 eV , -0.92918 eV 9 
-0.72457 eV , and -0.33582 eV for F, CI, Br , and /.respectively. 

» i i » * * 9 « * ■ 

Consider the case that the terminal bond is H- X corresponding to the angle of the 
atoms HCX where C is a carbon atom and X is a halogen atom. The factors c x and C x of 
Eq. (15.90) are 0,75 for all halogen atoms. For X = F, c' 2 is given by Eq. (15.69) with c 2 
of the participating carbon and F atoms given by Eq. (15.62) and Eq. (15.65), respectively. 
The factor C 2 of fluorine is one since it is the only halogen wherein the ionization energy is 
greater than that 13.605804 eV , the magnitude of the Coulombic energy between the electron 
and proton of H given by Eq. (1.243). For each of the other halogens, CT, Br, and /, cl is 
also given by Eq. (15.69) with c. of the participating carbon given by Eq. (15.62) and c, of 
the participating X atom given by c 2 = 0.91771 (Eq. (13.430)) since the X atom is energy 
matched. to the C2^ 3 HO. In these cases, C 2 is given by Eq. (15.65) for the corresponding 
atom X where C 2 matches the energy of the atom X to that of H . 
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Using the distance between the two atoms A and B of the general molecular 
group ACB when the total energy of the corresponding A - B MO is zero, the corresponding 
bond angle can be determined from the law of cosines: 

s t 2 + s 2 2 - 2^ 2 cosine 9 = s 3 2 (15 .96) 

With Sj = 2c' c ^ A , the internuclear distance of the C- A bond, s 2 = 2c' c _ B , the intemuclear 
distance of each C-B bond, and s y = 2c\_ B , the intemuclear distance of the two terminal 
atoms, the bond angle B M between the C - A and C - B bonds is given by 

(2*'c- J + (2c' C - J -2(2c« c ^X^^>sine^= (2*^ J (15.97) 



.(15.98) 



Consider the exemplary structure C i C a (O a )O b wherein C a is bound to C b , O o , and O y In 
the general case that the three bonds are coplanar and two of the angles are known, say 9 { and 
0 2 , then the third 0 3 can be determined geometrically: 

03 = 360-0,-0, (15.99) 
In the general case that two of the three coplanar bonds are equivalent and one of the angles is 
known, say Q x , then the second and third can be determined geometrically: 

(360- ft) 

g 2=^ = V 2 ' (15.100) 



ANGLES AND DISTANCES FOR AN MO THAT FORMS AN ISOSCELES 
TRIANGLE 

In the general case where the group comprises three A- B bonds having B as the central 
atom at the apex of a pyramidal structure formed by the three bonds with the A atoms at the 
base in the xy-plane. The C )v axis centered on B is defined as the vertical or z-axis, and any 
two A- B bonds form an isosceles triangle. Then, the angle of the bonds and the distances 
from and along the z-axis are determined from the geometrical relationships given by Eqs. 
(13.412-13.416): 

the distance d artg ^ B from the origin to the nucleus of a terminal B atom is given by 
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(15.101) 



ori£bt ' B 2sin60° 

the height along the z-axis from the origin to the A nucleus d is given by 



^ htlghi 



the angle d r of each A - B bond from the z-axis is given by 



g^tan^f ^^l (15.103) 



Consider the case where the central atom B is further bound to a fourth atom C and 
the B - C bond is along the z-axis. Then, the bond B tiABC given by Eq. (14.206) is 



DIHEDRAL ANGLE 

Consider the plane defined by a general ACA MO comprising a linear combination of;two 
C- A -bond MOs where C is the central atom. The dihedral angle 0 ZBCUCA between the 
ACA -plane and a line defined by a third bond with C , specifically that corresponding to a 
C - 2? -bond MO y is calculated from the bond angle and the distances between the A , 

5, and C atoms. The distance d x along the bisector of 0 ZAXA from C to the internuclear- 

distance line between A and A 9 2c % A _ A , is given by 

4 = 2c' c ^cos^^ . (15.105) 

where 2c % C _ A is the internuclear distance between A and C. The atoms A, A, and B 

define the base of a pyramid. Then, the pyramidal angle can be solved. from the 

internuclear distances between A and A, 2c \_ A , and between A and B, 2c ' AB , using the 
law of cosines (Eq. (15.98)); 



9 ZABA = C0S 



K 2c 'a-bX 2c 'a-b) 



(15.106) 



Then, the distance d t along the bisector of d, AM from B to the intermiclear-distance line 
2c is given by 
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4 = ^-b cos ^ 05.107) 

The lengths d x> d 2 , and 2c' c _ fl define a triangle wherein the angle between d x and the 

intemuclear distance between B and C 9 2c' C B9 is the dihedral angle 0^,^ that can be 
solved using the law of cosines (Eq. (1 5.98)): 



0 &CUCA =COS 



(15.108) 



SOLUTION OF GEOMETRICAL AND ENERGY PARAMETERS OF 
MAJOR FUNCTIONAL GROUPS AND CORRESPONDING ORGANIC 
MOLECULES 

The exemplary molecules given in the following sections were solved using the solutions of 
organic chemical functional groups as basis elements wherein the structures and energies 
where linearly added to achieve the molecular solutions. Each functional group can be treated 
as a building block to form any desired molecular solution from the corresponding linear 
combination. Each functional group element was solved using the atomic orbital and hybrid 

orbital spherical orbitsphere solutions bridged by molecular orbitals comprised of the H 2 - 

type prolate spheroidal solution given in the Nature of the Chemical Bond of Hydrogen-Type 
Molecules section. The energy of each MO was matched at the HO or AO by matching the 

hybridization and total energy of the MO to the AOs and HOs. The energy E (e.g. given 
by Eq. (15.58)) for a C2sp 3 HO and Eq.(15.59) for an 02p AO) was subtracted for each set 

■ 

of unpaired electrons created by bond breakage. 

The bond energy is not equal to the component energy of each bond as it exists in the 
molecule; although, they are close. The total energy of each group is its contribution to the 
total energy of the molecule as a whole. The determination of the bond energies for the 
creation of the separate parts must take into account the energy of the formation of any 
radicals and any redistribution of charge density within the pieces and the corresponding 
energy change with bond cleavage. Also, the vibrational energy in the transition state is 
dependent dn the other groups that are bound to a given functional group. This will effect the 
functional-group energy. But, because the variations in the energy based on the balance of the 
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molecular composition are typically of the order of a few hundreds of electron volts at most, 
they were neglected. 

The energy of each functional-group MO bonding to a given carbon HO is 
independently matched to the HO by subtracting the contribution to the change in the energy 
of the HO from the total MO energy given by the sum of the MO contributions and 

E(C 9 2sp*^= -14.63489 eV (Eq. (13.428)). The intercept angles are determined from Eqs. 

(15.71-15.78) using the final radius of the HO of each atom. The final carbon-atom radius is 
determined using Eqs. (15.32) wherein the sum of the energy contributions of each atom to all 
the MOs in which it participates in bonding is determined. This final radius is used in Eqs. 
(15.19) and (15.20) to calculate the final valence energy of the HO of each atom at the 
corresponding final radius. The radius of any bonding heteroatom that contributes to a MO is 
calculated in the same manner, and the energy of its outermost shell is matched to that of the 
MO by the hybridization factor between the carbon-HO energy and the energy of the 
heteroatomic shell. The donation of electron density to the AOs and HOs reduces the energy. 
The donation of the electron density to the MO's at each AO or HO is that which causes the 

resulting energy to be divided equally between the participating AOs or HOs to achieve 

• ■ > 

energy matching. 

The molecular solutions can be used to design synthetic pathways and predict product 

yields based on equilibrium constants calculated from the heats of formation. New stable 

■ . ■■ 

compositions of matter can be predicted as well as the structures of combinatorial chemistry 

reactions. Further important pharmaceutical applications include the ability to graphically or 

computationally render the structures of drugs that permit the identification of the biologically 

active parts of the molecules to be identified from the common spatial charge-density 

functions of a series of active molecules. Drugs can be designed according to geometrical 

: • ' 

parameters and bonding interactions with the data of the structure of the active site of the 

< 

drag. 

To calculate conformations, folding, and physical properties, the exact solutions of the 
charge distributions in any given molecule are used to calculate the fields, and from the fields, 
the interactions between groups of the same molecule or between groups on different 
molecules are calculated wherein the interactions are distance and relative orientation 

dependent. The fields and interactions can be determined using a finite-element-analysis 

.... 

approach of Maxwell's equations. 
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AROMATIC AND HETEROCYCLIC COMPOUNDS 

Aromatic and heterocyclic molecules comprise at least one of an aromatic or a cyclic 
conjugated alkene functional group. The latter was described in the Cyclic and Conjugated 
Alkenes section. The aromatic bond is uniquely stable and requires the sharing of the electrons 
of multiple H 2 -type MOs. The results of the derivation of the parameters of the benzene 

molecule given in the Benzene Molecule ( C 6 H 6 ) section can be generalized to any aromatic 
function group(s) of aromatic and heterocyclic compounds. 

C 6 J? 6 can be considered a linear combination of three ethylene molecules wherein a 

C-H bond of each CH 2 group of H 2 C = CH 2 is replaced by a C = C bond to form a six- 
member ring of carbon atoms. The solution of the ethylene molecule is given in the Ethylene 
Molecule (CH 2 CH 2 ) section. The radius r , (0.85252a 0 ) of the C2sp 2 shell of ethylene 

calculated from the Coulombic energy is given by Eq. (14.244). The Coulombic energy 

E a>utmb (pxhyiene'lsp 3 ^ (-15.95955 eV ) of the outeT electron of the C2sp 3 shell is given by Eq. 

(14.245). The energy E(^ emm ,2sp 3 ^ (-15.76868 e7)ofthe outer electron of the C2sp 3 shell 

is given by Eq. (14.246). E T (c = C,2^> 3 ) (-1.13380 eV) (Eq. (14.247), the energy change of 

each C2sp 3 shell with the formation of the C = C-bond MO is given by the difference between 
^(c^^Jlsp 3 ^ and E(C 9 2sj?y, C 6 H 6 can be solved.using the same principles as those used 

to solve ethylene wherein the 2s and 2p shells of each C hybridize to form a single 2sp 3 shell 

as an energy minimum, and the sharing of electrons between two C2sp 3 hybridized orbitals 
(HOs) to form a molecular orbital (MO) permits each participating hybridized orbital to 
decrease in radius and energy. Each 2sp 3 HO of each carbon atom initially has four impaired 
electrons. Thus, the 6 H atomic orbitals (AOs) of benzene contribute six electrons and the six 
sp 3 -hybridized carbon atoms contribute twenty-four electrons to form six C-H bonds and six 
C = C bonds. Each C-H bond has two paired electrons with one donated from the H AO 
and the other from the C2sp 3 HO. Each C = C bond comprises a linear combination of a factor 

* * * i * 

»• * ■ • 

of 0.75 of four paired electrons (three electrons) from two sets of two C2.sp 3 HOs of the 

participating carbon atoms. Each C-H and each C = C bond comprises a linear combination 
of one and two diatomic H 2 -type MOs developed in the Nature of the Chemical Bond of 

Hydrogen-Type Molecules and Molecular Ions section, respectively. 

Consider the case where three sets of C = C-bond MOs form bonds between the two 
carbon atoms of each molecule to form a six-member ring such that the six resulting bonds 
comprise eighteen paired electrons. Each bond comprises a linear combination of two MOs 

wherein each comprises two C2sp 3 UOs and 75% of a # 2 -type ellipsoidal MO divided between 
the C2sp 3 HOs: 

3 (2 C25p 3 + 0.75 H 2 MC>y> 3 (p = c\ ethylene - type - bond MO 
( * \ 

K-±6\C=C/- bond MO of benzene 



(15.142) 
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The linear combination of each # 2 -type ellipsoidal MO with each C2sp l HO further comprises 

an excess 25% charge-density contribution per bond from each C2sp 3 HO to the C = C-bond 
MO to achieve an energy minimum. Thus, the dimensional parameters of each bond C = C- 
bond are determined using Eqs. (15.42) and (15.1-15.5) in a form that are the same equations as 
those used to detennine the same parameters of the C = C-bond MO of ethylene (Eqs. (14.242- 
14.268)) while matching the boundary conditions of the structure of benzene. 

Hybridization with 25% electron donation to each C = C-bond gives rise to the 

,2sp> HO-shell Coulombic energy E Cadmb (c^,2^ 3 ) given by Eq. (14.245). To meet 
the equipotential condition of the union of the six C2sp 3 HOs, c 2 and C 2 of Eq. (15.42) for the 



bcittsnc 



3e 



aromatic C=C-bond MO is given by Eq. (15.62) as the ratio of 15.95955 eV 9 the magnitude of 
Ecou^fCtem*,,^ 3 ) (E* 04-245)), and 13.605804 eV 9 the magnitude of the Coulombic 
energy between the electron and proton of H (Eq. (1.243)): 

C 2 (benzeneC2sp 3 HO)= c 2 (benzeneC2sp* HO)= 13605804 eV = 0.85252 (15.143) 

15*95955 

. The, energies , of each C=C bond of benzene are also determined using the same 
equations as those of ethylene (Eqs. (14.251-14.253) and (14.319-14.333) with the parameters of 

benzene. Ethylene serves, as a basis element for the C=C bonding of benzene wherein each of 

3tf 

the six C=C bonds of benzene comprises (0.75)(4)= 3 electrons according to Eq. (15.142). 

- u 

The total energy of the bonds of the eighteen electrons of the C=C bonds of benzene, 
E t (&#6> c==c } is given by (6)(o.75) times E T+osc (c = c) (Eq. (14.492)), the total energy of 
the C=C -bond MO of benzene including the Doppler term, minus eighteen times E(c,2sp 3 ) 

3c 

(Eq. (14.146)), the initial energy of each C2^ HO of each C that forms the C=C bonds of 

bond order two. Thus, the total energy of the six C=C bonds of benzene with three electron 
per aromatic bond given by Eq. (14.493) is 

E T {cA>C=c) - (6X°-75K + «.,c(C = C)-(6X3>(C 2sp 3 ) 

= (6)(0.75X-66.05796eK)-18(-14.63489eF) (15.144) 
= -297.26081 eV- (-263.42798 eK) 
. = -33.83284 eV 

The results of benzene can be generalized to the class of aromatic and heterocyclic compounds. 
E hv of an aromatic bond is given by E T (ff 2 ) (Eqs. (1 1.212) and (14.486)), the maximum total 

* 

energy of each H 2 -type MO such that 

£17 1 . 



-31.63536831 eV K \ — + - 

Mc 2 2 V/* J 



(15.145) 



The factor of 0.75 corresponding to the three electrons per aromatic bond of bond order two 
given in the Benzene Molecule (C 6 H 6 ) section modifies Eqs. (15.52-15.56). Multiplication of 
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the total energy given by Eq. (15.55) by /, = 0.75 with the substitution of Eq. (15.145) gives the 
total energy of the aromatic bond: 



f E(basis energies) + E T (atom - atom, msp % ,AO^ 



E T (group)- f x 



-31.63536831 eV 



2fi4l 2- 



m 



mc 



2 



(15.146) 



The total bond energy of the aromatic group E D £™p) is the negative difference of the total 
energy of the group (Eq. (15.146)) and the total energy of the starting species given by the sum 
ofc A E m Jc 4 AO/HO)and c $ E MM (c 5 AOI HO): , 

f ( E(basis energies) + E T (atom - atom, msp* .AO) 



E D (?«w/>)- - 



-31.63536831 eV 



i 



2U 



\a C 2 e 2 



m 



mc 

e 



(15.147) 



Hc t E^(AO/HOyc s E^(c 5 AOJHO^ 

• * * i 

Since there are three electrons per aromatic bond, c A is three times the number of aromatic 
bonds. 

-. . Benzene can also be considered as comprising chemical bonds between six CH radicals 
wherein each radical: comprises a chemical bond between carbon and hydrogep atoms. The 
solution of the parameters of CH is given in the Hydrogen Carbide (CH) section. Those of the 
benzene are given in the Benzene Molecule (C 6 H 6 ) section. The energy components of V e , V p9 

T, V m> and E T are the same as those of the hydrogen carbide radical, except that 

E T (C = C,25p 3 )= -4.13379 eV (Eq. (14.247)) is subtracted from E T (CH) ofEq. (13.495)to 

match the energy of each C - H -bond MO to the decrease in the energy of the corresponding 

« 

C2sp 3 HO. In the corresponding generalization of the aromatic CH group, the geometrical 
parameters are determined using Eq. (15.42) and Eqs. (15.1-15.5) with 

E T (rtom-atom,msp*AOy=z -1.13379 eV . 

The total energy of the benzene C - H -bond MO, E T (C - H] 9 given by Eq. (14.467) 

is the sum of 0.5E T (c = C.lsp^ the energy change of each C2sp l shell per single bond due to 

the decrease in radiuis with the formation of the corresponding C=C -bond MO (Eq. (14.247)), 
and 4 E T (Cif), the a MO contribution given by Eq. (14:441). In the corresponding 

generalization of the aromatic CH group, the energy parameters are determined using Eqs. 

/ * -1.13379 eV _ _ 
(15.146-15.147) with = ! and E T [atom- atom,msp\AOy= . Thus, the energy 
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sum 



3* 

contribution to the single aromatic CH bond is one half that of the C=C double bond 

3c 

contribution. This matches the energies of the CH and C-C aromatic groups, conserves the 
electron number with the equivalent charge density as that of £ = 1 in Eqs. (15.18-15.21), and 

3c 

further gives a minimum energy for the molecule. Breakage of the aromatic C-C bonds to 
give CH groups creates unpaired electrons in these fragments that corresponds to c 3 = 1 in Eq. 

(1 5.56) with £ " given by Eq. (15.58). 

Each of the C - H bonds of benzene comprises two electrons according to Eq. (14.439). 

From the energy of each C - H bond, -E D ( l2 Ctf) (Eq. (14.477)), the total energy of the 

twelve electrons of the six C - H bonds of benzene, E T (C 6 H 6 ,C - H^ 9 given by Eq. (14.494) 
is 

E T (C 6 H 6 ,C-H)= (6)(-E d ^ (^CH^= 6(-3.90454 eV)= -23.42724 eV (15.148) 

The total bond dissociation energy of benzene, E D (c^H^ given by Eq. (14.495) is the negative 

( • le \ 

i of E T \C 6 H € ,C=C J (Eq. (14.493)) and E T {C 6 H 6> C - H) (Eq. (14.494)): 

■ = -((-33.83284 eV)± (-23.42724 eV^) (15.149) 
= 57.2601 eV 

Using the parameters given in Tables 15.214 and 15.216 in the general equatioiis (Eqs. (15.42), 
(15.1-15.5), and (15.146-15.147)) reproduces the results for benzene given in the Benzene 

Molecule (C 6 H 6 ) section as shown in Tables 15.214 and 15.216. 

The symbols of the functional groups of aromatics and hertocyclics are given in Table 
15.213. The geometrical (Eqs. (15.1-15.5) and (15.42)), intercept (Eqs. (I5.7i-15.78)), and 
energy (Eqs. (15.6-15.11), (15.17-15.56), and (15.146-15.147)) parameters of aromatics and 
hertocyclics are given in Tables 15.214, 15.215, and 15.216, respectively. The total energy of 
benzene given in Table 15.217 was calculated as the sum over the integer, multiple of each 

E D favup) of Table 15.216 corresponding to functional-group composition of the molecule. The 

bond angle parameters of benzene determined using Eqs. (15.79-15.108) are given in Table 
15.218. 

Table 15.213. The symbols of functional groups of aromatics and hertocyclics. 



Functional Group | Group Symbol ] 



CC (aromatic bond) C=C 
CH (aromatic) CH (i) 
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Table 15.214.. The geometrical bond parameters of aromatics and hertocyclics and 



Parameter 


3e 

c=c 

Group 


CH 
Group 




1.47348 


1.60061 


c> to 


1.31468 


1.03299 


Bond Length 2c' (□) 


1.39140 


1.09327 


Exp. Bond Length 

00 


1.399 
(benzene) 


1.101 

(benzene) 


b,c (a 0 ) 


0.66540 


1.22265 


e 


0.89223 


0.64537 



Table 15.216. The energy parameters (eV) of functional groups of aromatics and 



Parameters 

> 


it 

c=c 

Group 


CH 
Group 


fx. 


0.75 


1 




2 


1 




0 


0 




0 


0 




0.5 


0.75 

* 




0.85252 


1 




1 


1 




0.85252 


0.91771 




0 


1 


• 


3 


1 




0 


1 




0.5 


0.75 




0.85252 


1 


V t (eV) 


-101.12679 


-37.10024 


V p (eV) 


20.69825 


13.17125 


T(eV) 


34.31559 


11.58941 


V m XeV) 


-17.15779 


-5.79470 


EQoiho) (eV) 


0 


-14.63489 


^h^o W 1 ™) (eV) 


0 


-1.13379 
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E T {toiHo){eV) 


0 


-13.50110 


E T d, lM o)(eV) 


-63.27075 


-31.63539 


E T {atom - atom 9 msp 3 .AOj (eV) 


-2.26759 


-0.56690 


±j T {Mo)(ey) 


-65.53833 


-32.20226 


co HO 15 radls\ 


49.7272 


26.4826 




32.73133 


17.43132 


Id W 


-0.35806 


-0.26130 




0.19649 
[49] 


n or coo 

0.35532 
Eq. (13.458) 




-0.25982 


-0.08364 


E (eH 


0.14803 


0.14803 


E T (group) (e V) 


-49.54347 


-32.28590 




-14.63489 


-14.63489 


Km* 0*0 


0 


-13.59844 


£ D ^«mp) (eV) 


5.63881 


3.90454 
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1 Q 


0.79597 


0.79597 




|— 

ON 

• 


1= 

»— t 

o 


Final Total 
Energy 

C2sp 3 
(eV) 


fO 

oo 

DO 

CO 
IA 
1— < 
1 


CO 
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OO 

* 

CO 

»n 
»— < 
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o 
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vc 
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in 
ro 
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«n 
oo 

9 
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in 
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»n 

00 

» 
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oo 
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00 
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»— 1 




o" 








0.79232 








0.79232 




Atom 2 
Hybridization 
Designation 

(Table 
15 3 A) 


co 




•a 

1 B 

* 


-17.17218 


* 


Atoml 
Hybridization 
Designation 

(Table 
15.3 JV) 


to 




a 

1- 

SJ 


oo 

* 

f-H 
1 




2c' 

Terminal 


> ,/~\ 
1 Q 


v> 

00 

wo 
vo 




2c' 

Bond 2 


— s 

«° 


2.62936 




2c' 

Bond 1 


' — \ 


2.62936 




Atoms of 
Angle 


of 

81 


ZCCH 
(aromatic) 





120 [50-52] 

(benzene) 


120 [50-52] 
(benzene) 


o 


120.19 


119.91 












120.19 










vo 

CO 
00 

<r> 

00 

■ 

r— • 
1 










• 

o 




o~ 






Atoms of Angle 


zccc 

(aromatic) 


ZCCH 
(aromatic) 
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Software Program 

The present invention relates to a system of computing and rendering the nature of at 
least one specie selected from a group of diatomic molecules having at least one atom that is 
other than hydrogen, polyatomic molecules, molecular ions, polyatomic molecular ions, or 
molecular radicals, or any functional group therein, comprising physical, Maxwellian solutions 
of charge, mass, and current density functions of said specie, said system comprising: 
processing means for processing physical, Maxwellian equations representing charge, mass, 
and current density functions of said specie; and an output device in communication with the 
processing means for displaying said physical, Maxwellian solutions of charge, mass, and 
current density functions of said specie. 

In one embodiment, for example, the system comprises five components: (1) the 
graphical user interface (GUI); (2) the routine for parsing between an input chemical structure 
or name and taking the input and activating a routine to call up the parts of the molecule 
(functional groups), which are used for determining the energies and structure to be rendered; 
(3) the functional-group data base that has an organization of the theoretical solutions; (4) the 
rendering engine, which calculates and enables manipulations of the image, such as a three- 
dimensional model in response to commands, as well as responds to commands for data 

parameters corresponding to the image such as bond energies and charge distribution and 

* * * ■ * ■ ' •',,.< 

geometrical parameters; and (5) data transfer system for inputting numerical data into or out of 

■ • •:■!.* . • . • 

the computational components and storage components of the main system. The system 

further comprises spreadsheets with solutions of the bond parameters with output in any 

» • 

standard spreadsheet format. The system also comprises a data-handling program to, transfer 
data from the spreadsheets into the main program. 

The output may be, for example, at least one of graphical, simulation, text, and 
numerical data. The output may be the calculation of at least one of: (1) a bond distance 
between two atoms; (2) a bond angle between three of the atoms; (3) a bond energy between 
two atoms; (4) orbital intercept distances and angles; and (5) charge-density functions of 

• ■ * 

atomic, hybridized, and molecular orbitals, wherein the bond distance, bond angle, and bond 

energy are calculated from physical solutions of the charge, mass, and current density 

■ • • . . ■ • 

functions of atoms and atomic ions, which solutions are derived from Maxwell's equations 
using a constraint that a bound electron(s) does not radiate under acceleration. 

In other embodiments, the charge, current, energy, and geometrical parameters are 
output to be inputs to other programs that may be used in further applications. For example, 
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the data of heats of formation may be input to another program to be used to predict stability 
(existence of compounds) equilibrium constants and to predict synthetic pathways. That is, a 
novel composition of matter may be discovered by calculating at least one of a bond distance 
between two of the atoms, a bond angle between three of the atoms, and a bond energy 
between two of the atoms, orbital intercept distances and angles, charge-density functions of 
atomic, hybridized, and molecular orbitals, the bond distance, bond angle, and bond energy 
being calculated from physical solutions of the charge, mass, and current density functions of 
atoms and atomic ions, which solutions are derived from Maxwell's equations using a 
constraint that a bound electron(s) does not radiate under acceleration. The charge and current 
density functions may be used to predict the electric and magnetic fields of the species to 

determine other properties due to . the interaction of the fields between species. These fields 

■ . * • * i 

* . • * . * 

and the predictions of field interactions may be computed using Maxwell's equations. In one 
embodiment, finite-element analysis is used to predict or calculate the interaction and resulting 
properties, such as the freezing point, boiling point, density, viscosity, and refractive index. 
Furthermore, the output data can be used to give thennodynamic, spectroscopic, and other 

****** ► • * * 

properties, aid in drug design and other applications with or without direct visualization. 
Furthermore, the data can be input into other programs of the system, which calculate 
thermodynamic and other properties, or performs a simulation, such as a chemical reaction or 
molecular dynamics. 

* ■ i t ' * 

The output data may be used to predict a composition of matter comprising a plurality 
of atoms, the improvement comprising a novel property or use discovered by calculation of at 

■ • ' • i * * , . 

i . * ► i ■ 

least one of a bond distance between two of the atoms, a bond angle between three of the 
atoms, and a bond energy between two of the atoms, orbital intercept distances and angles, 
charge-density functions of atomic, hybridized, and molecular orbitals, the bond distance, bond 
angle, and bond energy being calculated from physical solutions of the charge, mass, and 

■ ♦ 

current density functions of atoms and atomic ions, which solutions are derived from 
Maxwell's equations using a constraint that a bound electron(s) does not radiate under 
acceleration. The novel property, for example, may be a new pharmaceutical use, or stability 
at room temperature of a novel arrangement of atoms or ions. 

In one embodiment, the output device of the system is a display that displays at least 

■ 

one of visual or graphical media. The display may be at least one of static or dynamic. At 
least one of vibration, rotation, and translation may be displayed. The displayed information 
may be used for at least one of modeling reactivity, predicting physical properties, and aiding 
in drug and material design. The output device may be a monitor, video projector, printer, or 

* » ■ ■ i 

■ * • * * • * 1 
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one-, two- or three-dimensional rendering device. The displayed information may be used to 

» * 

model other molecules and provides utility to anticipate their reactivity and physical properties. 
Additionally, data may be output and used in the same and additional applications as the 
rendered models and representations of the calculated physical solutions. The processing 
means of the system may be a general-purpose computer. The general-purpose computer may 
comprise a central processing unit (CPU), one or more specialized processors, system memory, 
a mass storage device such as a magnetic disk, an optical disk, or other storage device, an input 
means. The input means may comprise a serial port, USB port, microphone input, camera 
input, keyboard or mouse. The processing means comprises a special purpose computer or 
other hardware system. The system may comprise computer program products such as 
computer readable medium having embodied therein program code means. The computer 
readable media may be any available media which can be accessed by a general purpose or 
special purpose computer. The computer readable media may comprise, for example, at least 
one of RAM, ROM, EPROM, CD ROM, DVD or other optical disk storage, magnetic disk 
storage or other magnetic storage devices, or any other medium that can embody the desired 
program code means and which can be accessed by a general purpose or special purpose 

computer. The program code means may comprise executable instructions and data, which 

• > •• ■ . . . 

cause a general purpose computer or special purpose computer to perform a certain function of 

> , . *• < >. ■ * 

a group of functions. Commercial examples of suitable program language includes, for 

example, C++, C, JAVA, FORTRAN, Python and Assembly Languages, programmed with an 

• \ ■■• ....... . . • • 

algorithm based on the physical solutions, and the computer may be a PC, mainframe, 

supercomputer, or cluster of computers. Commercial examples of suitable programs include, 

■ 

for example, APIs like OpenGL, DirectX, FOX GUI toolkit, and Qt. This program may be 

* * » 

developed to run on at least one of operating systems like Windows XP, Windows 2000, 
Windows Vista, MAC OS, MAC OS X, Linux, Unix, Irix and other Unix-type operating 
systems. 

■ • 

Millsian software is designed to render 3-D models of molecules, molecular ions, 

... > 

molecular radicals, functional groups thereof, and related structure and property information 

• » 

and produce useful data output and application of the parameters of these species, wherein the 
nature of their bound electrons and chemical bonds are solved using Dr. Randell L. Mills' 

* * ■ 

Classical Quantum Mechanics theory described in the reference: R. L. Mills, "The Grand 
Unified Theory of Classical Quantum Mechanics", June 2006 Edition, Cadmus Professional 
Communications-Science Press Division, Ephrata, PA, ISBN 0963517171, Library of 
Congress Control Number 2005936834; posted at 
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http ://www.blackUghtpower.c»i^ookdownl^ which is incorporated by this reference 
in its entirety. Figure 76 provides a flow chart diagram, which is an example of a software 
system that can be utilized for this purpose, which example is not intended to limit the scope of 
the disclosed inventions. The main parts of this exemplary software system illustrated in 

4 

Figure 1 will now be further explained: 

Start: A user can start the program by running an executable program file. That might 
be done, for example, by double clicking the program icon on a Windows-based operating 
system, or typing the name of the executable file on the command line and pressing the 'Enter' 
key on a Linux or Unix operating system. The program initially starts by reading data files 
located in specific directories. The names of those directories and locations are fixed 
according to the type and format of the data files. 

Data files: There are two types of data files used in the software system: functional 
group data files and molecule data files. Functional group data files contain information about 
various functional groups. Functional groups are the basic bonding elements or units that each 
typically comprise an atom, or at least two atoms bound together as found within a molecule 
(e.g. -CI , C=C, C=0, CH3). Functional groups typically dictate or define properties and 
structure of the molecule. Similar functional groups in different molecules typically react in 
similar ways when subjected to a particular set of reaction conditions. Molecule data files 
contain information about molecules, molecular ions, and molecular radicals. These data files 
are processed according to their file formats. t 

File formats: Millsian software employs two kinds of file formats for storing 
information about the structure, energies and names of molecules and functional groups: raw- 
data format and hierarchical format. In a raw-data file, all information is stored as is, below 
the header describing the type of information. As shown in the Table below, for example, the 
names of the molecule or functional group are listed below the #NAMES header. The names 
and positions of the atoms are listed below the #ATOMS header, and so on. 
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if* t t, n ii in 

#NAMES 




CH3 




Alkane CH3 




#SMILES 




C- 




#ATOMS 




1 C 0 0 


0 


2 H 1.9775 -0.698 


0 


3 H -0.9888 


-0.698 1.712 


4 H -0.9888 


-0.698 -1.712 


#AO 




1 1 Is 0.171 2 




2 1 2sp3 0.864 -1 

• 




#BONDS 




1 12 1 1.649 




2 13 1 1.649 




3 14 1 1.649 




#BONDAXES 




11 1 90 


0 0 


#DATA 


* 


RCH bond angle: 109.44° 


• 


CH bond length: 2.097 a 


» 


HFORM: 12.492 





Most of the functional groups files and some simple molecule files are stored in raw 

4 

data format. The other file format used with the Millsian software system, the hierarchical file 
format, represents the information in a graph style, in which nodes are connected to other 
nodes through links. This format is designed to construct molecules by attaching different 
functional groups in a desired manner. Using this file format, the user can construct complex 
molecules built from solved functional groups. 

* 

For example, as shown in the Table below, the pentane data file includes a #GROUPS 
header below which is list all of the functional groups that form a part of the pentane molecule. 
Under the #GROUP_LINKS header is information about how these functional groups are 
connected to each other to construct pentane. 
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#NAMES 
Pentane 

#SMILES 

ccccc 

#GROUPS 

1 C- 

2 -C- 

3 -C- 

4 -C- 

5 C- 



#GROUP_LINKS 



1 


1 1 


2 1 


180 


2 


22 


3 1 


180 


3 


32 


41 


180 


4 


42 


5 1 


180 



Processing data files: As further shown in the flow diagram of Figure 1 for the 
exemplary Millsian software system, the program first processes the functional group data file 
and constructs the functional-group objects, which arc complete 3-D representations of the 
functional groups and their related information supplied in the corresponding file. These 
objects are then ready to be visually displayed through use of a molecule viewer. Next, the 
program processes the molecule data file. If the molecule date file is in raw-data format, then 

I 0 

the program makes a molecule object directly from it If the file is in hierarchical format, then 
the program calculates geometric parameters from listed functional groups. The program 
stores all functional-group objects and molecule objects using internal data structures. 

Visualization/ User Interactions: As shown in Figures 2 and 3, the molecule viewer 
displays the functional-group objects and molecule objects and provides basic interaction 
capabilities with the displayed objects, such as rotating, scaling, and moving the objects. The 
molecule viewer also provides other visualization options, such as viewing molecules in wire 
frame mode, viewing coordinate axes, and changing of the transparency and lighting. The 
user, for example, can also select parts of a molecule for visualization, like a nucleus, atomic 
orbital, molecular orbital, or bond axis. The viewer also includes a drop down information 
window, which provides, for example, related information about molecules, such as bond 

angles, component functional groups, and total heat of formation. In one embodiment, a user 

■ 

can create new molecules by joining functional groups. For this purpose, the user can select a 
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first functional group. Next, the user can select an open bond from the functional group where 
the user desires to attach another group. Next, the user can select another or the same 
functional group, followed by selecting an open bond from the other group. The user can join 
the two selected functional groups at the selected open bonds by clicking on 'Join Groups'. 
This method of joining functional groups at open bonds can be repeated to form the desired 
molecule. 

SMILES input The Simplified Molecular Input Line Entry Specification or SMILES is 
a specification for unambiguously describing the structure of chemical molecules using short 
ASCII strings. Through the user interface, a user can enter SMILES to construct molecules, 
provided the new molecule comprises functional groups that are in die database. Once a 
SMILES is entered, a parser reads in and breaks it down into component functional groups. 
The software system then attaches the component functional groups to create the new molecule 

... * * » t . 

* ... 

object, which can then be viewed using the molecule viewer. 

■ t ■ .>*•■• 

While the claimed invention has been described in detail and with reference to specific 

* * * * 

embodiments thereof, it will be apparent to one of ordinary skill in the art that various changes 
and modifications can be made to the claimed invention without departing from the spirit and 

■ 

scope thereof. 
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CLAIM: 

1 . A system of computing and rendering the nature of at least one specie selected from a 
group of diatomic molecules having at least one atom that is other than hydrogen, 
polyatomic molecules, molecular ions, polyatomic molecular ions, or molecular radicals, 
or any functional group therein, comprising physical, Maxwellian solutions of charge, 
mass, and current density functions of said specie, said system comprising: 

processing means for processing physical, Maxwellian equations representing charge, 
mass, and current density functions of said specie; and 

an output device in communication with the processing means for displaying said 
physical, Maxwellian solutions of charge, mass, and current density functions of said specie. 

♦ 

2. The system of claim 1 wherein the output device is a display that displays at least one of 

visual or graphical media, 

* * 

3. The system of claim 2 wherein the display is at least one of static or dynamic. 

4. The system of claim 3 wherein at least one of vibration and rotation is be displayed. 

5. The system of claim 4 wherein displayed information is used to model reactivity and 
physical properties. 

6. The system of claim 5, wherein the output device is a monitor, video projector, printer, or 
three-dimensional rendering device. 

7. The system of claim 6 wherein displayed information is used to model other species and 
provides utility to anticipate their reactivity and physical properties. 

8. The system of claim 7 wherein the processing means is a general purpose computer. 

9. The system of claim 8 wherein the general purpose computer comprises a central 
processing unit (CPU), one or more specialized processors, system memory, a mass 
storage device such as a magnetic disk, an optical disk, or other storage device, an input 
means. 
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10. The system of claim 9, wherein the input means comprises a serial port, usb port, 
microphone input, camera input, keyboard or mouse. 

1 1 . The system of claim 1 0 wherein the processing means comprises a special purpose 
computer or other hardware system. 

12. The system of claim 1 1 further comprising computer program products. 

13. The system of claim 12 further comprising computer readable medium having embodied 
therein program code means. 

i ■ » 

! 1 * . ' " I ■ . ■ . 

14. The system of claim 1 3 wherein the computer readable media is any available media 
which can be accessed by a general purpose or special purpose computer. 

* 

15. The system of claim 14 wherein the computer readable media comprises at least one of 
RAM, ROM, EPROM, CD ROM, DVD or other optical disk storage, magnetic disk 
storage or other magnetic storage devices, or any other medium which can embody the 
desired program code means and which can be accessed by a general purpose or special 
purpose computer. 

16. The system of claim 15 wherein the program code means comprises executable 
instructions and data which cause a general purpose computer or special purpose 
computer to perform a certain function of a group of functions. 

17. The system of claim 16 wherein the program code is Mathematica programmed with an 
algorithm based on the physical solutions, and the computer is a PC. 

18. The system of claim 17 wherein the algorithm is ParametricPlot3D[{2*Sqrt[l- 
z*z]*Cos[u],Sqrt[(l-z*z)]*Sin[u],z},{u,0,2* Pi},{z,-1,.9999}], and the rendering is 
viewed from different perspectives. 
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19. The system of claim 18 wherein the algorithms for viewing from different perspectives 
comprises Show[Out[l], ViewPoint->{x,y,z}] where x, y 5 and z are Cartesian 
coordinates. 

20. The system of claim 19 wherein the physical, Maxwellian solutions of the charge, mass, 
and current density functions of said specie comprises a solution of the classical wave 



21. The system of claim 20 wherein the boundary constraint of the wave equation solution is 
nonradiation according to Maxwell's equations. 



22. The system according to claim 1, wherein the force generalized constant k 1 of a H 2 -type 



where C, is the fraction of the H 2 -type ellipsoidal molecular orbital basis function of a 
chemical bond of the specie and C 2 is the factor that results in an equipotential energy match 
of the participating at least two molecular or atomic orbitals of the chemical bond. 

23. The system according to claim 1, wherein the distance from the origin of the MO to each 
focus c % is given by: 



equation V 2 -i-^ pQr,4A0- 
L v dt J 




ellipsoidal MO due to the equivalent of two point charges at the foci is given by: 




(15.1) 



4xe 0 




(15.2) 



the intemuclear distance is 




(15.3) 



the length of the semirainor axis of the prolate spheroidal MO ft = c is given by 



(15.4) 
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and, the eccentricity, e , is 

(15.5). 



e= — 
a 



24. The system according to claim 23, wherein a potential energy of the electrons in the 
central field of the nuclei at the foci is 



where ^ is the number of equivalent bonds of the MO for functional groups and in the case 
of independent MOs not in contact with the bonding atoms, the terms based on charge are 
multiplied by , the bond-order factor, which is 1 for a single bond, 4 for an independent 

double bond and 9 for an independent triplet bond. 

25. The system according to claim 1, wherein the potential energy of the two nuclei is 

V p =n i 05.7). 



26. The system according to claim 1, wherein the kinetic energy of the electrons is 



_ h 2 , ct+y/a 2 ~b 2 /icon 
T = w,c 1 c, 1 In , (15.8) 

where «, is the number of equivalent bonds of the MO for functional groups and in the case 
of independent MOs not in contact with the bonding atoms, the terms based on charge are 
multiplied by c m , the bond-order factor, which is 1 for a single bond, 4 for an independent 

double bond and 9 for an independent triplet bond. 

27. The system according to claim 1, wherein the energy, V m , of the magnetic force between 
the electrons is 
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where ^ is the number of equivalent bonds of the MO for functional groups and in the case 
of independent MOs not in contact with the bonding atoms, the terms based on charge are 
multiplied by , the bond-order factor, which is 1 for a single bond, 4 for an independent 

double bond and 9 for an independent triplet bond. 

28. The system according to claim 1 , wherein total energy of the H 2 -type prolate spheroidal 
MO, E T (h&o) , is given by the sum of the energy terms: 

E T [HM=Vt+T+V m +V p (15.10) 



%7teja 2 -b 2 



c,c 2 



I a) a-Ja 2 -b 2 



8«b*'L V a) a-c* J 



(15.11) 



where ^ is the number of equivalent bonds of the MO for functional groups and in the case 
of independent MOs not in contact with the bonding atoms, the terms based on charge are 
multiplied by , the bond-order factor, which is 1 for a single bond, 4 for an independent 

double bond and 9 for an independent triplet bond. 



29. The system according to claim 1, wherein the total energy E T [atom, msp 3 ) ( m is the 

integer of the valence shell) of the AO electrons and the hybridized shell is given by the 
sum of energies of successive ions of the atom over the n electrons comprising total 
electrons of the at least one AO shell according to the formula: 

E T (atom,msp 3 ) = -^IP a (15.12 
where IP m is the m th ionization energy (positive) of the atom. 



30. The system according to claim 29, wherein the radius of the hybridized shell is given 
by: 



. , . § Hgzfg (15.13). 
q-z-n ^SqEj. [atom, msp 3 j 
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3 1 . The system according to claim 29, wherein the Conlombic energy E CeuJomb {atom, msp* ) 
of the outer electron of the atom msp* shell is given by 



32. The system of claim 31, wherein in the case that during hybridization at least one of the 
spin-paired AO electrons is unpaired in the hybridized orbital (HO), the energy change for 
the promotion to the unpaired state is die magnetic energy E(magnetic) at the initial 

radius r of the AO electron: 

JIW0? 2B£L!aji. 05.15) 

then, the energy E(atom 9 msp 3 ) of the outer electron of the atom msp* shell is given by the 
sum of E Coulomb (atom,rnsp*) and E(magnetic): 

V 8 ^oV rn e r 

33. The system according to claim 29, wherein the total energy E T (moljatom,msp*) (m is 

the integer of the valence shell) of the HO electrons is given by the sum of energies of 
successive ions of the atom over the n electrons comprising total electrons of the at least 
one initial AO shell and the hybridization energy: 

n 

E T (molxitom,msp 3 ) = E[atom,msp 3 )-^IP m (15.17) 

m-2 

where lP m is the m th ionization energy (positive) of the atom and the sum of -IP X plus the 



hybridization energy is E[atom 9 msp*} . 



34. The system of claim 33, wherein the radius r, of the hybridized shell is given by: 
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V 2 (Z~q)-s(0.2S)) f ^ (15.18) 

^=z-rt J 8;rc 0 i? r [molMtom, msp* J 

where .y = 1,2,3 for a single, double, and triple bond, respectively. 

35. The systan of claim 34, wherein the Coulombic energy E Cotdomb {mol. atom, msp 3 ) of the 
outer electron of the atom msp 3 shell is given by 

E couiou* (molatom, msp 3 ) = (15.19). 

36. The system of claim 35, wherein in the case that during hybridization at least one of the 
spin-paired AO electrons is impaired in the hybridized orbital (HO), the energy change for 
the promotion to the unpaired state is the magnetic energy E(magnetic) at the initial 

radius r of the AO electron given by Eq. (1 5. 1 5). Then, the energy E {molatom, msp 3 ) of 

the outer electron of the atom msp 3 shell is given by the sum of E^,^ (mol.atom,msp 3 ) 
and E(magnetic) : 

E (moLatom, msp' ) = — — + ^f, (1 5 .20) 

E T {atom - atom, msp 1 ) , the energy change of each atom msp 3 shell with the formation of 
the atom-atom-bond MO is given by the difference between E (moLatom, msp 3 ) and 
E{atom,msp 3 ) : 

Ej [atom - atom, msp 3 ) = E{molxitom, msp 3 )-£ {atom, msp 3 ) (1 5.21). 

j 

37. The system of claim 3 1 , wherein E Coulomb {moLatom, msp 3 ) is one of: 

E Couhnb { C etye>* > 2sp 3 ) > E CouJomb [ C ethane > 2s P* ) > E Cou!omb [ C acetylene > 2s P* ) > 
^Coulomb {palkane » ) > 

£ co^6 [atom, msp 3 ) is one of 

E(mol.atom > msp 3 ) is one of £(C e(Meo ,,25p 3 ), Efc^.lsp 3 ) , 
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E(C^2sp>) E(C aikane9 2sp>); 

E{atom,msp 3 ) is one of and E(C>2sp 3 ) mdE(Cl,3sp 3 ); 
E T {atom~atam,msp 3 ) is one of E(C-C,2sp 3 ), E(C = C 9 2sp 3 ), and. 

E(C = C,2sp 3 )\ 

atommsp 3 is one of C2sp 3 9 Cl3sp 3 
E T {atom-atom(s x ),msp 3 ) is E T (C-C 9 2sp 3 ) and E r [atom- atom (s 2 ) 9 msp 3 ) is 

E T (c = C,2sp 3 ) y and 

In the case of the C2.srp 3 HO, the initial parameters (Eqs. (14.142-14.146)) are 



P -V ( 2 ~"> g2 - = ■ 10g2 c = 0.91771a 0 (15.22) 

^ ±£tej,(el48.25751eK) 8*e 0 (el48.25751 eV) 



-e 2 -e 1 



E r ,^(C,2sp 3 ) = — = = -14.82575 eV (15.23) 

Coi " om * 1 ' F > %xef w 8/K? 0 0.91771a 0 

Eimagnetic)=l!®&r= , =0,19086 «F (15.24) 

^(n) (0.84317*%) 

-e 2 In^h 2 



7W 



= -14.82575 eF +0.19086 (15.25) 
= -14.63489 eV 

wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to + 10%, if desired. 

38. The system of claim 37, wherein 

]T (Z-g)=10 (15.26) 

q=Z-n 

Equations (14.147) and (15.17) give 

^ (molMtom,msp l ) = E T (C e(toB ,2jp 3 ) = -151.61569 eV (15.27) 
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and using Eqs. (15,18-15.28), the final values of r cl ^ , E^^ (C2sp 3 ), and E(C2sp 3 ) , and 



the resulting E 7 



( B0 A 

\C-C 9 C2sp 3 



of the MO due to charge donation from the HO to the MO 



BO 



where C - C refers to the bond order of the carbon-carbon bond for different values of the 
parameter s are given in Table 1 5. L 



Table 15.1: 



MO 
Bond 
Order 
(BO) 


s 1 


s 2 


W ("o) 
Final 


^ OTi (C2^ 3 )(eV) 
Final 


E(C2sp 3 ) 
(eV) 


E T (c-C,C2sp 3 j 
(eV) 


I 


l 


0 


0.87495 


-15.55033 


-15.35946 


-0.72457 


n 


2 


0 


0.85252 


-15.95955 


-15.76868 


-1.13379 


m 


3 


0 


0.83008 


-16.39089 


-16.20002 


-1.56513 


IV 


4 


0 


0.80765 


-16.84619 


-16.65532 


-2.02043 



wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to + 10%, if desired. 

39. The system of claim 33, wherein a minimum-energy bond with the constraint that it must 
meet the energy matching condition for all MOs at all HOs or AOs, the energy 

E{moLatom y msp 3 ) of the outer electron of the atom msp 3 shell of each bonding atom is 
the average of E(molxitom 9 msp 3 ) for two different values of s : 



. x EimoLatom (s t ) , msp 3 ) + E imoljatom ( s 2 ) , msp 3 ) 
E [moLatom, msp J = — — — — 



(15.28). 



40. The system of claim 39, wherein in the case, E T {atom - atom, msp 3 ) , the energy change 

of each atom msp 3 shell with the formation of each atom-atom-bond MO, is average for 
two different values of s : 
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„ / A 3 n E T (atom-atom^^msp^^Eylatom-atomlsA^msp 3 ^ 

E T [atom - atom, msp 3 ) = — * i-ii — ~-L -± — L-L (15.29). 

2 

41 . The system of claim 39, wherein a first MO and its HOs comprising a linear combination 
of bond orders and a second MO that shares a HO with the first. In addition to the mutual 
HO, the second MO comprises another AO or HO having a single bond order or a mixed 
bond order, in order for the two MOs to be energy matched, the bond order of the second 
MO and its HOs or its HO and AO is a linear combination of the terms corresponding to 
the bond order of the mutual HO and the bond order of the independent HO or AO, and in 
general, E T [atom - atom, msp 3 ) , the energy change of each atom msp 3 shell with the 

formation of each atom-atom-bond MO, is a weighted linear sum for different values of s 
that matches the energy of the bonded MOs, HOs, and AOs: 

E T [atom - atom, msp 3 ) = £ c Sn E T [atom - atom (s n ) , msp 3 ) (1 5.3 0) 

where is the multiple of the BO of s n . The radius r^ 3 of the atom msp 3 shell of each 

bonding atom is given by the Coulombic energy using the initial energy E Couhmb [atom, msp 3 ) 

and E T [atom - atom, msp 3 ) , the energy change of each atom msp 3 shell with the formation 
of each atom-atom-bond MO: 

m rf o / 77 ( ^ 3\ . 1? ( T vvi" (15.31) 



Sft£ o a 0 J^c^mb (atom, msp 3 ) -f E T [atom - atom, msp 2 )) 



where E Coalomb (C2sp 3 ) = -14.825751 eV . 



42. The system of claim 41, wherein the Coulombic energy E CoM [molAtom,msp 3 } of the 
outer electron of the atom msp 3 shell is given by Eq. (15.19). 

43. The system of claim 42, wherein in the case that during hybridization, at least one of the 
spin-paired AO electrons is unpaired in the hybridized orbital (HO), the energy change for 
the promotion to the unpaired state is the magnetic energy E(magnetic) (Eq. (15.15)) at 

the initial radius r of the AO electron. Then, the energy E[molxitom,msp 3 ) of the outer 
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electron of the atom msp* shell is given by the sum of E Couloab {moLatom, msp 3 ) and 
E(magnetic) (Eq. (15.20)). 

♦ 

44. The system of claim 42, wherein, E 7 [atom - atom 9 msp 3 ), the energy change of each 
atom msp 1 shell with the formation of the atom-atom-bond MO is given by the difference 
between E(molxitom,msp 2 ) and E[atom,msp 3 ) given by Eq. (15.21). 



45. The system of claim 44, wherein using the equation (15.23) for E^^ (C, 2sp* ) in 

equation (1 5 .3 1), the single bond order energies given by Eqs. (15.1 8-1 5.27) and shown in 
Table 15.1, and the linear combination energies (Eqs. (15.28-15.30)), the parameters of 
linear combinations of bond orders and linear combinations of mixed bond orders are 
given in Table 15.2, wherein the calculated and measured values and constants recited in 
the equations herein can be adjusted, for example, up to + 10%, if desired. 

46. The system of claim 29, wherein the radius r^^, of the C2.jp 3 HO of a carbon atom of a 

given specie is calculated using Eq. (14.5 14) by considering £ [M0 9 2sp 3 ) , the total 
energy donation to each bond with which it participates in bonding. 



47. The system of claim 1, wherein equation for the radius is given by 

= -V 

W 8^o [E CM (C,2sp*)+Y, ^ {MO, 2sp> )) 

e 2 



(15.32). 



8^r 0 (el4.825751 eK + £|^ (M0 9 2sp 3 )\) 



48. The system of 29, wherein the Coulombic energy E C(mlomb {molxitom 9 msp 3 } of the outer 
electron of the atom msp 3 shell is given by Eq. (1 5.19). 
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49. The system of claim 29, in the case that during hybridization, at least one of the spin- 
paired AO electrons is unpaired in the hybridized orbital (HO), the energy change for the 
promotion to the unpaired state is the magnetic energy E{magnetic) (Eq. (15.15)) at the 

- 

initial radius r of the AO electron. Then, the energy E [mo\jatom> msp* ) of the outer 

electron of the atom msp 3 shell is given by the sum of E Coutofab [molatom,msp* ) and 
E{magnetic) (Eq. (15.20)). 



50. The system of claim 49, wherein for the C2sp 3 HO of each methyl group of an alkane 

contributes -0.92918 eV (Eq. (14.513)) to the corresponding single C-C bond; the 

corresponding C2sp* HO radius is given by Eq. (14.514). The C2sp 3 HO of each 

methylene group of C n H 2ri+2 contributes -0.92918 eV to each of the two corresponding 

G - C bond MOs^ wherein the calculated and measured values and constants recited in 
the equations herein can be adjusted, for example, up to + 10%, if desired. 

51 . The system of claim 50, wherein the radius (Eq. (15.32)), the Coulombic energy (Eq. 
(15.19)), and the energy (Eq. (15.20)) of each alkane methylene group are 

2 



r _ ~e 



(Ecou** (C, V )+ Z £ w (methylene C - C,2sp 3 )) 



e 2 



lm<, (el4.825751 eF+e0.92918 eV+e0.929UeV) 
= 0.81549^, 



(15.33) 



E{C nethvl 2sp*) = -r^- -+ 2w - r = -16.49325 eV (15.35) 

\ V ) ^(O^C^) ^(Q.84317^) 3 ^ ' 

wherein the calculated and measured values and constants recited in the equations 
herein can be adjusted, for example, up to + 10%, if desired. 

52. The system of claim 5 1 , wherein in the determination of the parameters of functional 
groups, heteroatoms bonding to C2sp 3 HOs to form MOs are energy matched to the 
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Clsp 3 HOs, the radius and the energy parameters of a bonding heteroatom are given by 
the same equations as those for C2sp 3 HOs. 

53. The system of claim 52, wherein using Eqs. (15.15), (15.19-15.20), (15.24), and (15.32) 
in a generalized fashion, the final values of the radius of the HO or AO, r AtomHOAO , 

E Coulomb (^oLatom 7 msp 3 ) 9 and E(C mol 2sp 3 ) are calculated using J^E^MO^sp 3 ), 
the total energy donation to each bond with which an atom participates in bonding 

corresponding to the values of E T - C, C2sp 3 j of the MO due to charge donation 

from the AO or HO to the MO given in Tables 1 5. 1 and 1 5.2 and the final values of 
r AumMOAO , E Co ^{molMtom,msp 3 ),and E(C mol 2sp 3 ) calculated using the values of 

£ r ^C^C,C2sp 3 j given in Tables 15.1 and 15.2 are shown in Tables 15.3Aand 15.3B 

in the specification, wherein the calculated and measured values and constants recited in 
the equations herein can be adjusted, for example, up to ± 10%, if desired. 

54. The system of claim 53, wherein the energy of the MO is matched to each of the 
participating outermost atomic or hybridized orbitals of the bonding atoms wherein the 
energy match includes the energy contribution due to the AO or HO's donation of charge 
to the MO. 

55. The system of claim 54, wherein the force constant k % (Eq. (15.1)) is used to determine 
the ellipsoidal parameter c ' (Eq. (1 5.2)) of the each i/ 2 ~type-ellipsoidal-MO in terms of 
the central force of the foci, c % is substituted into the energy equation (from Eq. (15.1 1))) 
which is set equal to n, times the total energy of H 2 where n, is the number of equivalent 
bonds of the MO and the energy of H 2> -31.63536831 eV , Eq. (11.212) is the minimum 

energy possible for a prolate spheroidal MO, wherein the calculated and measured values 
and constants recited in the equations herein can be adjusted, for example, up to + 10%, if 
desired. 
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56. The system of claim 55, wherein the energy equation and the relationship between the 
axes, the dimensions of the MO are solved, the energy equation has the semimajor axis a 
as it only parameter, the solution of the semimajor axis a allows for the solution of the 
other axes of each prolate spheroid and eccentricity of each MO (Eqs. (15.3-15.5)), and 
the parameter solutions then allow for the component and total energies of the MO to be 
determined. 



57. The system of claim 1, wherein the total energy, E T (h^mo) , is given by the sum of the 
energy terms (Eqs. (15.6-15.1 1)) plus E T {AO I HO) : 



E T (H lM o) = V e +T+V m +V p +E T (AOIHO) (15.36) 

a + yla 2 -b 2 



Ej (HiM0) = 



\ a) a->Ja 2 - 



+ E T (AO/HO) 



= — ^- c^-^lln— "I +EJAOIHO) 

(15.37) 

where is the number of equivalent bonds of the MO, c, is the fraction of the H 2 -type 
ellipsoidal MO basis function of a chemical bond of the group, c 2 is the factor that results in 
an equipotential energy match of the participating at least two atomic orbitals of each 
chemical bond, and E T (AO J HO) is the total energy comprising the difference of the energy 

E(AO/HO) of at least one atomic or hybrid orbital to which the MO is energy matched and 

any energy component AE HtM0 (AO/HO) due to the AO or HQ's charge donation to the MO. 



58. The system of claim 57, wherein 

E T (AOfHO) = E(AO/HO)-AE H ^ 0 (AO/HO) (15.38) 



59. The system of claim 58, wherein as specific examples, E T (AO I HO) is one from the 
group of 

E T {AO I HO) = E(02p shell) = -EQonlzation; O) = -1 3.6 1 8 1 e V ; 
E T (AO I HO) = E(N2p shell) = -EQomzation; N) = -14.53414 eV ; 
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'Coulomb 



AO I HO) = E(c,2sp 3 ) = -14.63489 eV ; 

(Cl,3sp 3 ) = -14.60295 eV ; 
AO I HO) = EQonization; C) + E(iomzation; C+) ; 
AOIHO) = E (C alam ,2sp 3 ) = -l535946eV; 

AO/HO^+Efc^lsp^-Efc^lsp 3 ); 

AO I HO) = E(C, 2sp 3 ) - 2E T (C = C, 2sp 3 ) = -14.63489 eV -(-2.26758 eV) ; 

= ^(C^, V)~ ^(C^.,2^ 3 )- £(C^ ew ,2^ 3 ) = 16.20002 ; 
= E(C,2sp 3 ) -2E T (C a C,2sp 3 ) = -14.63489 eV -(-3.13026 eF) ; 
= E (C 6emtne ,2sp 3 ) - ^(C^, 2^) ; 

= E{c,2sp 3 )-E r (C =C,2sp 3 ) = -14.63489 eK-(-l. 13379 eK), and 



E T 


(AO /HO) 


E T 


(AO /HO) 


E T 


(AO /HO) 


E T 


(AO/HO) 


E T 


(AO /HO) 


E T 


(AO /HO) 


E T 


(AO /HO) 




(AO /HO) 




[AO/HO) 


E T \ 


[AO! HO) 




[AO/HO) 



values and constants recited in the equations herein can be adjusted, for example, up to + 
10%, if desired. 



> 60. The system of claim 1, wherein to solve the bond parameters and energies, 



c' = a 



o 



m/2C x C z a \2C X C 2 



(Eq. (15.2)) is substituted into E T (h&o) to give 



E T [h 2 uo) = - 



Sxs n yja 2 -b 2 



c x c 2 



. a+Va 2 -6 2 , 
In , -1 



+ E r (AO/HO) 



r 



%7t€ Q C 



c x c 2 



UE T {AOIHO) 



*0 



a j a-c 



rip* 



2QQ 



<?,C 2 



-1 



2C,C 2 



+ Ej.(AO/HO) 



(15.39) 



wherein the total energy is set equal to E(basis energies) which in the most general case is 
given by the sum of a first integer «, times the total energy of H 2 minus a second integer 
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times the total energy of H , minus a third integer ^ times the valence energy of E{AO) 

(e.g. E(N) = -14.53414 eV) where the first integer can be 1,2,3..., and each of the second 
and third integers can be 0,1,2,3.... 

61. The system of claim 60, wherein 

E(basis energies) = q (-3 1 .6353683 1 eV) - ^ (-13.605804 eV)- n,E(AO) 
(15.40) 

in the case that the MO bonds two atoms other than hydrogen, E(basis energies) is n, times 
the total energy of H 2 where «, is the number of equivalent bonds of the MO and the energy 
of H 2 y -31 .6353683 1 eV 9 Eq. (1 1 .212) is the minimum energy possible for a prolate 
spheroidal MO: 

E(basis energies) = i\ (-31.63536831 eV) (15.41) 

Ep (h 2 mo) , is set equal to E(basis energies) , and the semimajor axis a is solved. 

62. The system of claims 23 or 60, wherein the semimajor axis a is solved from the equation 
of the form: 



+ E T {AO I HO) = E{basis energies) 
(15.42) 

The distance from the origin of the H 2 -type-ellipsoidal-MO to each focus c 1 , the intemuclear 
distance 2c 1 , and the length of the semiminor axis of the prolate spheroidal # 2 -type MO 
b = c are solved from the semimajor axis a using Eqs. (15.2-15.4). 

63. The system of claim 62, wherein the component energies are given by Eqs. (15.6-15.9) 
and (15.39). 

* 

64. The system of claim 57, wherein the total energy of the MO of the functional group, 
E T (mo) , is the sum of the total energy of the components comprising the energy 
contribution of the MO formed between the participating atoms and 




89 



WO 2007/053486 



PCT/US2006/042140 



E T [atom - atom, msp l AO) , the change in the energy of the AOs or HOs upon forming 



the bond. 



65. The system of claim 64, wherein from Eqs. (15.39-15.40), E T (mo) is 

E T (mo) = E(basis energies) + E T {atom - atom, msp* AO) (1 5.43) 



66. The system of claim 65, wherein during bond formation, the electrons undergo a reentrant 
oscillatory orbit with vibration of the nuclei, and the corresponding energy E osc is the sum 



of the Doppler, E D , and average vibrational kinetic energies, E 



Kvtb • 



E CS c^n x {E D + E Kvib )^n l 



E i 



hv\\ 2 

mc 



(15.44) 



where «, is the number of equivalent bonds of the MO , k is the spring constant of the 
equivalent harmonic oscillator, and /i is the reduced mass. 

67. The system of claim 66, wherein the angular frequency of the reentrant oscillation in the 
transition state corresponding to E D is determined by the force between the central field 
and the electrons in the transition state. 



68. The system of claim 67,wherein the force and its derivative are given by 



f(R)=-c t 



'BO 



4xe 0 R' 



(15.45) 



and 



4xe 0 R : 



(15.46) 



such that the angular frequency of the oscillation in the transition state is given by 




C C e' 
4xe 0 R 3 



(15.47) 



where R is the semimajor axis a or the semiminor axis b depending on the eccentricity of 
the bond that is most representative of the oscillation in the transition state, c B0 is the bond- 
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order factor which is 1 for a single bond and when the MO comprises r\ equivalent single 
bonds as in the case of functional groups, c B0 is 4 for an independent double bond as in the 
case of the C0 2 and N0 2 molecules and 9 for an independent triplet bond, C lo is the 
fraction of the H 2 -type ellipsoidal MO basis function of the oscillatory transition state of a 
chemical bond of the group, and C 2o is the factor that results in an equipotential energy match 
of the participating at least two atomic orbitals of the transition state of the chemical bond. 
Typically, C l0 = C, and C 2o - C 2 , the kinetic energy, E K , corresponding to E D is given by 
Planck's equation for functional groups: 



69. The system of claim 68, wherein the Doppler energy of the electrons of the reentrant orbit 




is 




(15.49) 



given by the sum of E D and E KvA is 




\ 



E. 



{group) = n i (E D +E Xvi6 ) = n l E^ 



(15.50) 



Efa, of a group having r\ bonds is given by E T {mo)/^ such that 




is given by the sum of £ 7 M (Eq. (15.42))and E^ (Eq. (15.51)): 
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8 ^ lKI 

1 2C A 



V a J 



In 





1 aa 0 


2C,C 2 






2C,C 2 



1 



^(AO/ HO) + E T ( atom - atom, msp* A O) 













m.c 2 



1 . [ft 



= (£(&<m* energies) + E T [atom - atom, w^p 3 .^o)) 1 + 1 — * — ^ 



2^ 


'OA,,* 2 






m,c 2 

V 



+ 77,— ft 

^2 



(15.52). 



70. The system of claim 69, wherein the total energy of the functional group E T (group) is the 
sum of the total energy of the components comprising the energy contribution of the MO 
formed between the participating atoms, E(basis energies) , the change in the energy of 

the AOs or HOs upon forming the bond ( E T [atom - atom, msp*AO) ), the energy of 

oscillation in the transition state, and the change in magnetic energy with bond formation, 



E mag- 



71. The system of claim 70, wherein from Eq. (15.52), the total energy of the group E T (croup) 
is 



Ey (proup) = 



[E(basis energies) + Ej. [atom - atom, msp y -AOfj 1 + ^ 





\c l0 c 2 y 








m e c 2 



\ 



mag 



(15.53). 
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72. The system of claim 71, wherein the change in magnetic energy E mag which arises due to 

the formation of unpaired electrons in the corresponding fragments relative to the bonded 
group is given by 



"** 3 my * r 3 



(15.54) 



where r 3 is the radius of the atom that reacts to form the bond and % is the number of 



electron pairs. 
f 



Ey (Group) = 



^E(basis energies) + E T {atom - atom, msp^AO^ 



1+ 



\ 



C C e d 



2^ 



m. 



mc 

c 



3 J 

r 



(15.55). 

73 . The system of claim 72, wherein the total bond energy of the group E D (Group) is the 
negative difference of the total energy of the group (Eq. (15.55)) and the total energy of 
the starting species given by the sum of c k E mm (c 4 AO I HO) and c 5 E Jnflfal (c 5 AO/HO) : 



E D (Group) = 



{Eipasis energies) + E T (atom-atom,msp 2 ^AOf) 
+r h E Mb + c fi^-(c A E Mtud {AO I HO) + c 5 E Mael (c 5 AO I HO)) 

















TO, 





(15.56). 



74. The system of claim 70, wherein in the case of organic molecules, the atoms of the 
functional groups are energy matched to the C2sp 3 HO such that 

E(AO/HO) --14.63489 eV (15.57) 

10 For examples of E mQg from previous sections: 
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E m (C2sp 3 ) = c 3 l&A = c % JVA = c 3 0.14803 eV (15.58) 

{02p) = c 3 = c 3 ^1 = c 3 0.1 1441 eV (1 5.59) 

£ (N2p) = c 3 ^£ = c 3 3 = <r 3 0.14185 eV (15.60). 

' 3 r* 3 (0.93084^7, ) 3 



75. The system of claim 74, wherein in the general case of the solution of an organic 

functional group, the geometric bond parameters are solved from the semimajor axis and 
the relationships between the parameters by first using Eq. (15.42) to arrive at a , the 
remaining parameters are determined using Eqs. (15.1-15.5), the energies are given by 
Eqs. (15.52-15.59), and to meet the equipotential condition for the union of the # 2 -type- 

ellipsoidal-MO and the HO or AO of the atom of a functional group, the factor c 2 of a 

H 2 -type ellipsoidal MO in principal Eqs. (15.42) and (1 5.52) may given by 



(i) one: 

c 2 =l (15.61) 

(ii) the ratio that is less than one of 13.605804 eV , the magnitude of the Coulombic 
energy between the electron and proton of H given by Eq. (1.243), and the magnitude of the 

Coulombic energy of the participating AO or HO of the atom, E Coahmb {MO.atom, msp 3 ) 
given by Eqs. (15.19) and (15.31-15.32). For \E Cmtlonb (MO.atom 9 msp 3 )\>l3. 605&04 eV : 



e 2 



8ra? 0 q 0 _ 13.605804 eV 



fas*/" . _ . - 

0 A-B AorBsp 



E Cmkmt (MOMtom,msp>)\ 



(15.62) 



For 



E Coulomb {MOxaom, msp* )| < 1 3.605804 e V : 



e 2 



^Coulomb 



2 e 2 13.605804 eV 

&ae 0 a 0 



(15.63) 
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(iii) the ratio that is less than one of 13-605804 eV , the magnitude of the Coulombic 
energy between the electron and proton of H given by Eq. (1-243), and the magnitude of the 
valence energy, Eiyalence) , of the participating AO or HO of the atom where Eiyalence) is 

the ionization energy or E[MOxitom 9 msp 3 ) given by Eqs. (15.20) and (15.31-15.32). For 
\E(valence)\ >13.605804 eV : 



fay, ".MSB 04 eV 
e \E(yalence)\ 

For \E(valence)\< 13. 605804 eV : 
e 2 



m ^ m faience)] 

e 2 13.605804 eV 1 



(iv) the ratio that is less than one of the magnitude of the Coulombic energy of the 
participating AO or HO of a first atom, E Coulomb (MO.atom,msp*) given by Eqs. (15.19) and 

(15.31-15.32), and the magnitude of the valence energy, Eiyalence), of the participating AO 
or HO of a second atom to which the first is energy matched where Eiyalence) is the 

ionization energy or E^MOxttom 9 msp 3 ) given by Eqs. (15.20) and (15.31-15.32). For 

• * ■ 

\ E cauhmb [MO. atom, msp 3 )| > Eiyalence) : 

\E(yalenc4 ^ 

|^c~/«»i (MOxitom, msp 3 )| 

For \E Cmb>nb (MO.atom,msp 3 )\<E(valence) : 

\B Cmtmh (MQjitom t msp^ ^ ^ 
2 \E{vdlence)\ 

(v) the ratio that is less than one of the magnitude of the valenoe-Ievel energies, 
E n iyalence) , of the AO or HO of the nth participating atom of two that are energy matched 
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where Eiyahnce) is the ionization energy or E [MOjitom, msp 3 j given by Eqs. (15.20) and 

(15.31-15.32): 

= ^valence) 
E^iyalence) 

(vi) the factor that is the ratio of the hybridization factor c 2 (l) of the valence AO or 

HO of a first atom and the hybridization factor c 2 (2) of the valence AO or HO of a second 

atom to which the first is energy matched where c 2 (») is given by Eqs. (15.62-15.68); 

alternatively c 2 is the hybridization factor c 2 (1) of the valence AOs or HOs a first pair of 

atoms and the hybridization factor c 2 (2) of the valence AO or HO a third atom or second 
pair to which the first two are energy matched: 

c 2 =^ (15.69) 
c 2 {2) 



(vii) the factor that is the product of the hybridization factor c 2 (l) of the valence AO 
or HO of a first atom and the hybridization factor c 2 (2) of the valence AO or HO of a second 
atom to which the first is energy matched where c 2 (n) is given by Eqs. (15.62-15.69); 
alternatively c 2 is the hybridization factor c 2 (l) of the valence AOs or HOs a first pair of 

atoms and the hybridization factor c 2 (2) of the valence AO or HO a third atom or second 

pair to which the first two are energy matched: 
^=^(1)^(2) (15.70) 

The hybridization factor c 2 corresponds to the force constant k (Eqs. (1 1 .65) and (13.58)). 

In the case that the valence or Coulombic energy of the AO or HO is less than 13.605804 eV , 
the magnitude of the Coulombic energy between the electron and proton of H given by Eq. 
(1.243), then C 2 corresponding to k' (Eq. (15.1)) is given by Eqs. (15.62-15.70). 



76. The system of claim 75, wherein specific examples of the factors c 2 and C 2 of a # 2 -type 
ellipsoidal MO of Eq. (15.51) given in following sections are 
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C 4ayH OM F) = £^c,(C2^O) = 5^|£L (0 . 9]77l) = „. 77M:; 

c, (c^ho to aj-g^ (cv») - S«| = ..,,3.7= 
c;(eV«»..i)- 5 ^^(cV«>)- ^^g (Mi77i)-a«»7. 

/„; 1Hrt -14.53414 <?K ■ 

{H " 1N) ' l&kj " -15.35-M6.K - °- 94627 ; 

c, (« „ »,) = . . 0.93383 . 

' E[C,2sp 3 ) -15.56407 eV 

^ /« ' rr\ •£('S'>3p) -10.36001 erf 

c 2 (O to S3sp 3 to C2sp 3 HO)=^Qc 2 (C2sp 3 HO) 

-13.61806 eV 



-10.36001 eV 
= 1.20632 



(0.91771); 



„ /„- 3 \ Ecouionifasp 3 ) -11.57099 eV noenAC 
cJS3sp ) = ^ = =0.85045: 

n y) E{H) -13.60580 eV 

4 

Q(CV» ,o *V ) -^gL (<V ) - E^f^W . 0.66,5, , 
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C 2 (S3sp 3 to O to C2sp i H0) = y^y^ c 2 (C2sp 3 HO) 

-1 1.52126 eV 



-13.61806 eV 
= 0.77641 



(0.91771) ; 



c 2 (O to N2p to C2sp 3 HO) = Cj (cisp'HO) 

-13.61806 eV 



-14.53414 eV 
= 0.85987 



(0.91771); 



~ x c 2 (C2sp 3 HO to N) 0 91140 

c,(JV2/> to 02n) = -24 =Hii^ =] 06727 • 

2K P P) c 2 (C2sp 3 HOtoO) 0.85395 ' ' 

C 2 (bemeneC2sp 3 HO) = c 2 (benzemC2 S p 3 H6) = l3605S04eV = 085252 . 
c 2 (arylC2sp 3 HO to O) = E ^°) c 2 (ary!C2sp 3 HO) 



E(c,2sp 3 ) 

-13.61806 eV 

~ -14.63489 eV 
= 0.79329 



(0.85252) 



C 2 



(H to anlineA^ *<& -14.53414eF 
V ' E(c,2sp 3 ) -15.76868 eV uyzi/i ' 

c 2 (aryWsjtHO to N) = ?} N \, c 2 (arylC2sp 3 HO) 

EiC 9 2sp j ' 

-14.63489 eK v ' 
= 0.84665 

r xt\ E ( N ) -14.53414 eV 
c 2 (H to anhne A r )= > v ' < = = 0.92171 : 

' E(C,2sp 3 ) -15.76868 eV 
c 2 (arylC2sp 3 HO to N)= ffi^V ? ("tfC^iTO) 



£(C 3 2jp 3 ) 

_ -14.53414 eV 

~ -14.63489 eV 
= 0.84665 



(0.85252) , and 



C 



2 



foi , t*. r>n izrA\ E( S > 3 P) -10.36001 eV 
(S3p to aryl-type C2sp 3 HO 1= / = = 0. 

V ' E(C,2sp*) -15.76868 eV 



65700. 
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wherein the calculated and measured values and constants recited in the equations herein 
can be adjusted, for example, up to ± 10%, if desired. 



77. The system of claim 1 , wherein the parameters of the point of intersection of each H 2 - 
type ellipsoidal MO and the A -atom AO are determined from the polar equation of the 
ellipse: 

1+e 



r = r n 



l + ecostf* 



(15.71). 



78. The system of claim 77, wherein the radius of the A shell is r A , and the polar radial 
coordinate of the ellipse and the radius of the A shell are equal at the point of intersection 
suchthat 



r A =(a-c'} 



a 



(15.72). 



1 + -COS0' 

a 



79. The system of claim 78, wherein the polar angle 6 ' at the intersection point is given by 



0 f = cos 



-i 



a 

7 



1 + 



( a - c >\ flL_i 



(15.73). 



80. The system of claim 79, wherein the angle 0 AAO the radial vector of the A AO makes 
with the internuclear axis is 



(15.74). 



8 1 . The system of claim 80, wherein the distance from the point of intersection of the orbitals 
to the internuclear axis is the same for both component orbitals such that the angle 
cot = 6 HiMO between the internuclear axis and the point of intersection of each H 2 -1ype 
ellipsoidal MO with the A radial vector obeys the following relationship: 



r A S™&AAO= bsine HMO 



(15.75) 



such that 
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^■^ ikgfrg (15 .76). 



82. The system of claim 81, wherein the distance d H ^ 0 along the internuclear axis from the 
origin of H 2 -type ellipsoidal MO to the point of intersection of the orbitals is given by 
= acosd HiUO (15.77). 



83. The system of claim 82, wherein the distance d A A0 along the internuclear axis from the 

origin of the A atom to the point of intersection of the orbitals is given by 
d AA0 = c'-d U2U0 (15.78). 

* 

84. The system of claim 1, wherein in ACS MO comprising a linear combination of C-A- 
bond and C - B -bond MOs where C is the general central atom and a bond is possible 
between the A and B atoms ofthe C-A and C-JB bonds, the /ACB bond angle is 
determined by the condition that the total energy ofthe H 2 -type ellipsoidal MO between 
the terminal A and B atoms is zero. 



85. The system of claim 84, wherein the force constant k ' of a H 2 -type ellipsoidal MO due 
to the equivalent of two point charges of at the foci is given by: 

C C 2e 2 
k'= ^> Ze (15.79) 
4xe 0 

where C, is the fraction of the H 2 -typc ellipsoidal MO basis function of a chemical bond of 
the specie which is 0.75 (Eq. (13.59)) for a terminal A-H {A is H or other atom) and 1 
otherwise and C 2 is the factor that results in an equipotential energy match ofthe 

participating at least two atomic orbitals of the chemical bond and is equal to the 
corresponding factor of Eqs. (15.42) and (15.52). 

86. The system of claim 85, wherein the distance from the origin ofthe MO to each focus c 1 
of the A-B ellipsoidal MO is given by: 



a ^m e e 2 2qC 2 a V2C ] C 2 ( 15 - 80 )' 
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87. The system of claim 85, wherein the intemuclear distance is 



; 88. The system of claim 85, wherein the length of the semiminor axis of the prolate 
spheroidal A - B MO b = c is given by Eq. (15.4). 

89. The system of claim 85, wherein the component energies and the total energy, E T (h&o) , 
of the A-B bond are given by the energy equations (Eqs. (11.207-11.212), (11.213- 

1 1 .21 7), and (1 1 .239)) of H 2 except that the terms based on charge are multiplied by c w , 
the bond-order factor which is 1 for a single bond and when the MO comprises r\ 
equivalent single bonds as in the case of functional groups. c B0 is 4 for an independent 
double bond as in the case of the C0 2 and N0 2 molecules. 

90. The system of claim 85, wherein the kinetic energy term is multiplied by which is 1 
for a single bond, 2 for a double bond, and 9/2 for a triple bond, the electron energy terms 
are multiplied by c x , the fraction of the H 2 -type ellipsoidal MO basis function of a 
terminal chemical bond which is 0.75 (Eq. (13.233)) for a terminal A-H (A is H or 
other atom) and 1 otherwise. 

91. The system of claim 85, wherein the electron energy terms are further multiplied by c\ , 
the hybridization or energy-matching factor that results in an equipotential energy match 
of the participating at least two atomic oibitals of each terminal bond. 

92. The system of claim 84, wherein when A-B comprises atoms other than H> 

E T [atom - atom.msp^ .AO) , the energy component due to the AO or HO's charge 

donation to the terminal MO, is added to the other energy terms to give E T (h&o) : 



E T (ffjMO) = 



8xe 0 c 



2 




1 +E T (atom-atom 9 msp 3 .AO 



(15.82). 



101 



WO 2007/053486 



PCT/US2006/042140 



93. The system of claim 84, the radiation reaction force in the case of the vibration of A - B 
in the transition state corresponds to the Doppler energy, E D , given by Eq. (1 LI 81) that 

is dependent on the motion of the electrons and the nuclei. 

94. The system of claim 93, wherein the total energy that includes the radiation reaction of the 
A-B MO is given by the sum of Ej. (h&o) (Eq. (15.82)) and given Eqs. (1 1.2 13- 
11.220), (11.231-11.236), and (11.239-1 1.240). 



95. The system of claim 94, wherein the total energy E T (A - B) of the A-B MO including 
the Doppler term is 



E T (A-B) = 



r 5 
—e 



v 



8^ 0 c' 





f 




cd 


^BO U BO 









a-c' 



+ E T (atom-atom,msp 3 .A(?) 



V 
J 




1 c&e 2 


c *o<? 




8ns 0 (a+c'f 


M 



(15.83) 

where C lo is the fraction of the 7/ 2 -type ellipsoidal MO basis function of the oscillatory 
transition state of the A - B bond which is 0.75 (Eq. (13.233)) in the case of H bonding to a 
central atom and 1 otherwise, C 2o is the factor that results in an equipotential energy match of 
the participating at least two atomic orbitals of the transition state of the chemical bond, and 

H = 171 i m i i s the reduced mass of the nuclei given by Eq. (11.1 54). 



96. The system of claim 95, wherein to match the boundary condition that the total energy of 
the A - B ellipsoidal MO is zero, E T (A-B) given by Eq. (15.83) is set equal to zero and 
substitution of Eq. (15.81) into Eq. (15.83) gives 
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0 = 



\ 



c x c 2 



I aa 0 

^BO c BO \ Ui I 1 

V 2C,C 2 



atom - atom 



,msp 2 .AO} 



J 



l+ \ 



'BO 



1W 



m. 



'50 



glg2 g2 

8^ 0 a + 



\3 



aa 0 



V 2C,C 2 



(15.84). 



97. The system of claim 96, wherein the vibrational energy-term of Eq. (15.84) is determined 

by the forces between the central field and the electrons and those between the nuclei 

i 

(Eqs. (11.141-11.145)). 

98. The system of claim 96, wherein the electron-central-field force and its derivative are 
given by 



f(0) = -c B o J ? S2L T (15-85) 



and 



/'(«) = 2c M ^fL (15.86). 



4ffe 0 cr 



99. The system of claim 97, wherein the nuclear repulsion force and its derivative are given 
by 

2 



: /(a+c') = 



%jte 0 {a+c*y 



(15.87) 



and 



4^(a + c f ) : 



(15.88) 



such that the angular frequency of the oscillation is given by 



0) = 




(15.89). 
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1 00. The system of claim 99, wherein since both terms of E m =E D + E^ are small due to 

the large values of a and c\ an approximation of Eq. (15.84) which is evaluated to 
determine the bond angles of functional groups is given by 



0 = 



r 



2C,C 2 [ 



cA 2-^lln 



a+ 1 s - 



2C,C 2 



-1 



a — 



2qc 2 



J 



E T {atom - atom, msp 3 .A O) 



2^ 


1 c,e 2 






m e c 2 



+ 



1*1 



\3 



2Cf 2 



(15.90). 



101. The system of claim 100, wherein from the energy relationship given by Eq. (15.90) 
and the relationship between the axes given by Eqs. (1 5.2-1 5.5), the dimensions of the 
A-B MO can be solved. 

102. The system of claim 101, wherein Eq. (15.90) is solved by the reiterative technique 
using a computer. 

103. The system of claim 100, wherein a factor c 2 of a given atom in the determination of 
c\ for calculating the zero of the total A- B bond energy is given by Eqs. (15.6245.65). 



104. The system of claim 103, wherein in the case of a if -if terminal bond of an alkyl 
or alkenyl group, c\ is the ratio of c 2 of Eq. (1 5.62) for the H-H bond which is one and 

c 2 of the carbon of the corresponding C-H bond: 

1 13.605804 eV 



c 2 = 



c 2 (C2sp>) E CM (C-HC2sp>) 



(15.91) 



wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to + 10%, if desired. 
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105. The system of claim 100, wherein in the case of the determination of the bond angle 
of the ACH MO comprising a linear combination of C - A -bond and C-H -bond MOs 
where A and C are general, C is the central atom, and c 2 for an atom is given by Eqs. 

(15.62-15.70), c\ of the A -H terminal bond is the ratio of c 2 of the A atom for the 

A - H terminal bond and c 2 of the C atom of the corresponding C - H bond: 

d = c 2 (A(A-H)msp 3 ) 
2 c 2 (C(C-H)(™p 3 ) 

» 

106. The system of claim 100, wherein in the case of the determination of the bond angle 
of the COH MO of an alcohol comprising a linear combination of C - O -bond and 
O-H -bond MOs where C, £?, and H are carbon, oxygen, and hydrogen, respectively, 
c 2 of the C - H terminal bond is 0.91771 since the oxygen and hydrogen atoms are at 

the Conlomb potential of a proton and an electron (Eqs. (1.236) and (10.162), 

respectively) that is energy matched to the C2sp 3 HO, wherein the calculated and 

measured values and constants recited in the equations herein can be adjusted, for 
example, up to + 10%, if desired. 

* 

107. The system of claim 1 00, wherein in the determination of the hybridization factor 

of Eq. (15.90) from Eqs. (15.62-15.70), the Coulombic energy, E Couiomb (MO.atom,msp 3 ) , 

or the energy, E [MOxitom, msp* ) , the radius r A B A ^ of the A or B AO or HO of the 

heteroatom of the A - B terminal bond MO such as the C2sp z HO of a terminal C-C 
bond is calculated using Eq. (15.32) by considering ^E Tmi [MO,2sp 3 }> the total energy 
donation to each bond with which it participates in bonding as it forms the terminal bond. 

1 08 . The system of claim 1 07, wherein the Coulombic energy E CouUmb {MO.atom> msp 2 ) of 
the outer electron of the atom msp 3 shell is given by Eq. (15.19). 



109. The system of claim 108, wherein in the case that during hybridization, at least one of 
the spin-paired AO electrons is unpaired in the hybridized orbital (HO), the energy change 
for the promotion to the unpaired state is the magnetic energy E(magnetic) (Eq. (15.15)) 
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at the initial radius r ofthe AO electron, and the energy E(MO*am,mrf) oftheouter 

electron ofthe atom msp* shell is given by the sum of E Cmdomb [MO.atom, msp*) and 
E(magnetic) (Eq. (15.20)). 

110. The system of claim 99, wherein in the specific case of the terminal bonding of two 
carbon atoms, the factor of each carbon given by Eq. (15.62) is determined using the 

Coulombic energy E^^ (C-C C2sp 3 ) ofthe outer electron ofthe Clsp 3 shell given 
by Eq. (15.19) with the radius of each C2sp 3 HO of the terminal C-C bond 

calculated using Eq. (15.32) by considering ( MO, 2sp % ) , the total energy donation 

to each bond with which it participates in bonding as it fonns the terminal bond including 
the contribution ofthe methylene energy, 0.92918 eV (Eq. (14.513)), corresponding to 
the terminal C-C bond. The corresponding E T ^atom-atom > msp 3 ,AO)mEq. (15.90) is 

E T [C-C C2sp 3 ) = -1 .85836 eV , wherein the calculated and measured values and 

constants recited in the equations herein can be adjusted, for example, up to ± 10%, if 
desired. 

111. The system of claim 1 00, wherein in the case that the terminal atoms arc carbon or 
other heteroatoms, the terminal bond comprises a linear combination ofthe HOs or AOs; 
thus, c' 2 is the average of the hybridization factors ofthe participating atoms 
corresponding to the normalized linear sum: 

c\ =^(c' 2 (atom \)+c' 7 (atom 2)) (15.93) 
In the exemplary cases of C-C, 0-0, and N-N where C is carbon: 



2 2 



( 2 2 ^ 

e e 



e 2 e 2 



I 

2 



^^^A-AAAOIHO ^ ne tf A-A kiAOIHO ) 

( 13.605804 eV 13.605804 eV 



(15.94) 



\ E co„ h »> (A-A^AO/HO) (A-AA.AO/HO) 
In the exemplary cases of C-N , C-O, and C-S, 



J 
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1 13.605804 eV 



+c 2 (C toB) 



(15.95) 



where C is carbon and c 2 (C to B) is the hybridization factor of Eqs. (15.52) and (15.84) that 
matches the energy of the atom B to that of the atom C in the group. 

112. . The system of claim 97, wherein the corresponding E T {atom - atom, msp 3 AO} term 

in Eq. (15.90) depends on the hybridization and bond order of the terminal atoms in the 
molecule, but typical values matching those used in the determination of the bond 
energies (Eq. (15.56)) are 

E T (C - O C2sp\02p) = -1.44915 eV ; E, (C - O C2sp 3 .02p) = -1 .65376 eV ; 



E T (C -N C2sp 3 .N2p) = -1.44915 eV ; E 7 (C-S C2sp s S2p) = -0.72457 eV ; 
E r (0-0 02p.02p) = -1.44915 eV ; E, (O- O 02p.02p) = -1.65376 eV ; 
E T {N-NN2pmp)= -1.44915 eV; E J (N-0 N2p.02p) = -1AA9\5 eV; 
E T {F-FF2pJ?2p) = -\M9\5eV; E T (Cl-Cl CBpJCBp) = -0.92918 eV; 
E T {Br- Br Br4p.Br4p) = -0.92918 eV ; E T {I -I I5p.I5p) = -0.36229 eV ; 
E T {C-F C2sp } J?2p) = -l.Z5t36eV;E r {C-Cl C2sp\Cl3p) = -0.9291 8 eV; 
E T {C-Br C2sp\BrAp) = -0.72457 eV ; E T (C-7 C2sp\l5p)= -0.36228 eV , and 



JE r (O-C/O2/).a3/>) = -0.92918eF . . + . . 1+ . . 

v ' , wherem the calculated and measured values and 



constants recited in the equations herein can be adjusted, for example, up to + 10%, if desired. 

113. The system of claim 100, wherein in the case that the terminal bond is X- X where 
X is a halogen atom, c, is one, and c[ is the average (Eq. (15.93)) of the hybridization 
factors of the participating halogen atoms given by Eqs. (15.62-15.63) where 



E Coulomb (MO.atom, msp*) = 13.605804 eV for X = / . 
1 14. The system of claim 1 13, wherein the factor C x of Eq. (15.90) is one for all halogen 



E couiomb (MOMtom 9 msp*) is determined using Eq. (15.32) and 



atoms. 
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115. The system of claim 1 13, wherein the factor C 2 of fluorine is one since it is the only 

» 

halogen wherein the ionization energy is greater than that 13.605804 eV , the magnitude 
of the Coulombic energy between the dectron and proton of H given by Eq. (1 .243), 
wherein the calculated and measured values and constants recited in the equations herein 
can be adjusted, for example, up to + 10%, if desired 

1 16. The system of claim 113, wherein for each of the halogens, Cl,Br 9 and / , C 2 is 
the hybridization factor of Eq. (15.52) given by Eq. (15.70) with c 2 (l) being that of the 
halogen given by Eq. (1 5.68) that matches the valence energy of X (E x (valence) ) to that 
of the C2sp 3 HO (E 2 (valence) = -14.63489 eV,Eq. (15.25)) and to the hybridization of 

Clsp 1 HO (c 2 (2) = 0.91771 , Eq. (13.430)), wherein the calculated and measured values 

and constants recited in the equations herein can be adjusted, for example, up to + 10%, if 
desired. 



117. The system of claim 116, wherein E T [atom - atom, msp*AO) of Eq. (15 .90) is the 

maximum for the participating atoms which is -1.44915 eV , -0.92918 eV , 
-0.92918 eV , and -0.33582 eV for F , CI , Br , and / , respectively, wherein the 
calculated and measured values and constants recited in the equations herein can be 
adjusted, for example, up to + 10%, if desired. 

1 1 8. The system of claim 100, wherein in the case that the terminal bond is C - X where 
C is a carbon atom and X is a halogen atom, the factors c, and Q of Eq. (15.90) are 

one for all halogen atoms. 

119. The system of claim 118, wherein for X = F , c 2 is the average (Eq. (1 5.95)) of the 
hybridization factors of the participating carbon and F atoms where c 2 for carbon is 
given by Eq. (15.62) and c 2 for fluorine matched to carbon is given by Eq. (15.70) with 
c 2 (1) for the fluorine atom given by Eq. (15.68) that matches the valence energy of F 
(E l (valence) = -17.42282 eV) to that of the CZsp 3 HO {E 2 (vahnce) = -14.63489 eV, 
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Eq. (15.25)) and to the hybridization of C2sp> HO (c 2 (2) = 0.91771, Eq. (13.430)), 

wherein the calculated and measured values and constants recited in the equations herein 
can be adjusted, for example, up to + 1 0%, if desired. 

120. The system of claim 119, wherein the factor Cj of fluorine is one since it is the only 

halogen wherein the ionization energy is greater than that 13.605804 eV , the magnitude 
of the Coulombic energy between the electron and proton of H given by Eq. (L243), 
wherein the calculated and measured values and constants recited in the equations herein 
can be adjusted, for example, up to ± 10%, if desired. 

121 . The system of claim 1 19, wherein for each of the other halogens, CI , Br , and / , c\ 

is the hybridization fector of the participating carbon atom since the halogen atom is 
energy matched to the carbon atom. 

122. The system of claim 121, wherein the C 2 of the terminal-atom bond matches that 
used to determine the energies of the corresponding C - X -bond MO. 

123. Hie system of claim 122, wherein C 2 is the hybridization factor of Eq. (15.52) given 
by Eq. (15.70) with c 2 (l) for the halogen atom given by Eq. (1 5.68) that matches the 
valence energy of X (E^ {valence) ) to that of the C2sp 3 HO 

(E 2 (valence) = -14.63489 eV , Eq. (15.25)) and to the hybridization of C2sp y HO 

(Cj (2) = 0.91771 , Eq. (13.430)), wherein the calculated and measured values and 

constants recited in the equations herein can be adjusted, for example, up to + 10%, if 
desired. 

124. The system of claim 123, wherein E T (atom - atom, msp*AO) of Eq. (1 5.90) is the 

maximum for the participating atoms which is -1.85836 eV , -0.92918 eV , 
-0.72457 eV , and -0.33582 eV for F , CI , Br , and / , respectively, wherein the 
calculated and measured values and constants recited in the equations herein can be 
adjusted, for example, up to + 10%, if desired. 
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125. The system of claim 11 9, wherein in the case that the terminal bond is H-X 
corresponding to the angle of the atoms HCX where C is a carbon atom and if is a 
halogen atom, the factors c x and C, of Eq. (15.90) are 0.75 for all halogen atoms. 



126. The system of claim 1 1 % wherein for X = F , c 2 is given by Eq. (1 5.69) with c 2 of 
the participating carbon and F atoms given by Eq. (15.62) and Eq. (15.65), respectively. 

127. The system of claim 126, wherein the factor C 2 of fluorine is one since it is the only 

halogen wherein the ionization energy is greater than that 13.605804 eV , the magnitude 
of the Coulombic energy between the electron and proton of H given by Eq. (1.243), 

« * 

wherein the calculated and measured values and constants recited in the equations herein 
can be adjusted, for example, up to ± 10%, if desired. 

128. The system of claim 127, wherein for each of the other halogens, CI, Br 9 and / , c 2 
is also given by Eq. (15.69) with c 2 of the participating carbon given by Eq. (15.62) and 
c 2 of the participating X atom given by c 2 = 0.91771 (Eq. (13.430)) since the X atom is 
energy matched to the C2sp 3 HO. 



129. The system of claim 128, wherein C 2 is given by Eq. (15.65) for the corresponding 
atom X where C 2 matches the energy of the atom X to that of H . 



130. The system of claim 100, wherein the distance between the two atoms A and B of 
the general molecular group ACB when the total energy of the corresponding A - B MO 
is zero, the corresponding bond angle can be determined from the law of cosines: 
s } 2 -f s 2 - 2^ 2 cosine^ = s 3 2 (15.96). 



131. The system of claim 1 30, wherein with s x = 2c' c _ A , the internuclear distance of the 
C - A bond, s 2 = 2c' c _ B , the internuclear distance of each C-B bond, and = 2c } A _ B , 
the internuclear distance of the two terminal atoms, the bond angle 0^ CB between the 
C -A and C-B bonds is given by 
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(2 C ' c _,) 2 +(2c' c . J ,) J -2(2 C ' C -,)(2c' c . fl )cosine^ = (2c',. i) ) 2 (15.97) 



GjACB = COS 



f (2c' c .,) a + (2c' c , J? f-(2c',. i> ) 2 ^ 



2(2 C ^)(2c' c . B ) 



(15.98). 



132. He system of claim 131, wherein the structure C 0 C a (O a Y) b wherein C„ is bound to 
C b , 0 o , and O b , the three bonds are coplanar and two of the angles are known, say 0, 
and 0 2 , then the third 0 3 can be determined geometrically: 

i 

03=360-0,-02 (15.99) 

133. The system of claim 131, wherein in the general case that two of the three coplanar 
bonds are equivalent and one of the angles is known, say 0, , then the second and third can 
be determined geometrically: 

*-«-G«tal (1 5,oo). 

134. The system of claim 1 , wherein in the general case where the group comprises three 
A - B bonds having B as the central atom at the apex of a pyramidal structure formed by 
the three bonds with the A atoms at the base in the xy-plane. 

135. The system of claim 134, wherein the C 3v axis centered on B is defined as the 

vertical or z-axis, and any two A-B bonds form an isosceles triangle, and the angle of 
the bonds and the distances from and along the z-axis are determined from the 
geometrical relationships given by Eqs. (13.412-13.416): 
the distance d origin _ B from the origin to the nucleus of a terminal B atom is given by 

d ori . R = 2c *- s (15.101) 
arigftt ~ B 2 sin 60° 

the height along the z-axis from the origin to the A nucleus d heigftt is given by 

4** = >/( 2c '^) 2 ~( rf ^) 2 > 0 5. 102) 
the angle 0 V of each A - B bond from the z-axis is given by 
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9. = tan' 1 



(d. 



\ 



orfgin-B 



(15.103). 



136. The system of claim 1 35, wherein in the case where the central atom B is further 
bound to a fourth atom C and the B - C bond is along the z-axis. Then, the bond 
given by Eq. (14.206) is 

0*bc =180-0,. (15.104). 



137. The system of claim 1 34, wherein in the plane defined by a general ACA MO 
comprising a linear combination of two C-A -bond MOs where C is the central atom, 
the dihedral angle 0^ CMC4 between the ACA -plane and a line defined by a third bond 
with C , specifically that corresponding to a C - B -bond MO, is calculated from the bond 
angle Q^ca m & the distances between the A , B , and C atoms. 



138. The system of claim 137, wherein the distance d A along the bisector of 0^ from 
C to the internuclear-distance line between A and A , 2c \_ A , is given by 

4 = 2c' c ^cos^^- (15.105) 
where 2c * C _ A is the internuclear distance between A and C . 



1 39. The system of claim 138, wherein the atoms A , A , and B define the base of a 
pyramid. 

140. The system of claim 139, wherein the pyramidal angle can be solved from the 
internuclear distances between A and A , 2c * A _ A , and between A and B , 2c \_ B , using 
the law of cosines (Eq. (15.98)): 

f (2c^) 2 + (2 C y j ) 2 -(2cV,) : 
2{2c> A _ B ){2c> A _ B ) 



OzABA = COS"' 



(15.106). 
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141. The system of claim 140, wherein the distance d 7 along the bisector of 9 ZABA from B 
to the internuclear-distance line 2c , is given by 

d 2 = 2c\_ B cos^^ (15,107). 



142. The system of claim 141, wherein the lengths d 19 d 2 , and 2c * M define a triangle 
wherein the angle between and the intermiclear distance between B and C , 2c X C _ B , is 
the dihedral angle Qtnc IACA that can be solved using the law of cosines (Eq. (15.98)): 

(d? + (2c> c _ B ) 2 -dl) 



&ZBC/ACA ~ C0S 



I 24(2c' B _ c ) J 



(15.108). 



143. The system of claim 1, wherein the specie are solved using the solutions of organic 
chemical functional groups as basis elements wherein the structures and energies where 
linearly added to achieve the molecular solutions, each functional group can be treated as 
a building block to form any desired molecular solution from the corresponding linear 
combination, each functional group element was solved using the atomic orbital and 
hybrid orbital spherical orbitsphere solutions bridged by molecular orbitals comprised of 
the H 2 -type prolate spheroidal solution given in the Nature of the Chemical Bond of 

Hydrogen-Type Molecules section, the energy of each MO was matched at the HO or AO 
by matching the hybridization and total energy of the MO to the AOs and HOs, the energy 

E mag (e.g. given by Eq. (15.58)) for a C2sp* HO and Eq.(15.59) for an 02p AO) was 
subtracted for each set of unpaired electrons created by bond breakage. 

1 44. They system of claim 143 , wherein the bond energy is not equal to the component 
energy of each bond as it exists in the specie, although, they are close. 

■ 

145. The system of claim 143, wherein the total energy of each group is its contribution to 
the total energy of the specie as a whole. , 

146. The system of claim 143, wherein the determination of the bond energies for the 
creation of the separate parts must take into account the energy of the formation of any 
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radicals and any redistribution of charge density within the pieces and the corresponding 
energy change with bond cleavage. 

1 47. The system of claim 143, wherein the vibrational energy in the transition state is 
dependent on the other groups that are bound to a given functional group, which will 
effect the functional-group energy, however because the variations in the energy based on 
the balance of the molecular composition are typically of the order of a few hundreds of 
electron volts at most, they are neglected. 

148. The system of claim 143, wherein the energy of each functional-group MO bonding to 
a given carbon HO is independently matched to the HO by subtracting the contribution to 
the change in the energy of the HO from the total MO energy given by the sum of the MO 

contributions and E(C,2sp 3 ) = -14.63489 eV (Eq. (13.428)). 

149. Tie system of claim 143, wherein the intercept angles are determined from Eqs. 
(15.71-15.78) using the final radius of the HO of each atom. 

1 50. The system of claim 143, wherein a final carbon-atom radius is determined using Eqs. 
(15.32) wherein the sum of the energy contributions of each atom to all the MOs in which 
it participates in bonding is determined. 

■ 

151. The system of claim 150, wherein the final radius is used in Eqs. (15.19) and (15.20) 
to calculate the final valence energy of the HO of each atom at the corresponding final 
radius. 

1 52. The system of claim 151, wherein the radius of any bonding heteroatom that 
contributes to a MO is calculated in the same manner, and the energy of its outermost 
shell is matched to that of the MO by the hybridization factor between the carbon-HO 
energy and the energy of the heteroatomic shell. 

153. The system of claim 150, wherein the donation of electron density to the AOs and 
HOs reduces the energy. 



114 



WO 2007/053486 



PCT/US2006/042140 



1 54. The system of claim 143, wherein the donation of the electron density to the MO's at 
each AO or HO is that which causes the resulting energy to be divided equally between 
the participating AOs or HOs to achieve energy matching. 

1 55. The system of claim 1, wherein the molecular solutions are used to design synthetic 
pathways and predict product yields based on equilibrium constants calculated from the 
heats of formation. 

1 56. The system of claim 1, wherein the new stable compositions of matter are predicted as 
well as the structures of combinatorial chemistry reactions. 

157. The system of claim 1 , wherein pharmaceutical applications include the ability to 
graphically or computationally render the structures of drugs that permit the identification 
of the biologically active parts of the specie to be identified from the common spatial 
charge-density functions of a series of active species. 

158. The system of claim 1 , wherein novel drugs are designed according to geometrical 
parameters and bonding interactions with the data of the structure of the active site of the 
drug. 

1 59. The system of claim 1 , wherein to calculate conformations, folding, and physical 
properties, the exact solutions of the charge distributions in any given specie are used to 
calculate the fields, and from the fields, the interactions between groups of the same 
specie or between groups on different species are calculated wherein the interactions are 
distance and relative orientation dependent. 

160. The system of claim 159, the fields and interactions can be determined using a finite- 
element-analysis approach of Maxwell's equations. 

1 61 . The system of claim 1 , wherein in the case where three sets of C = C -bond MOs 
form bonds between the two carbon atoms of each molecule to form a six-member ring 
such that the six resulting bonds comprise eighteen paired electrons, and each bond 
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comprises a linear combination of two MOs wherein each comprises two C2sp 3 HOs and 
75% of a H- -type ellipsoidal MO divided between the C2sp 3 HOs: 



3(2 C2sp 3 +0.75H 2 MOy* 3(e=c\ ethylene -type- bond M(?\ 



(15.142). 



bond MO of benzene 



1 62. The system of claim 1 62, wherein the linear combination of each H 2 -type ellipsoidal 
MO with each C2sp 3 HO further comprises an excess 25% charge-density contribution 
per bond from each C2sp 2 HO to the C = C -bond MO to achieve an energy minimum. 

. . . . . 

163. The system of claim 1 63, wherein the dimensional parameters of each bond C = C - 

■ 

bond are determined using Eqs. (15.42) and (15.1-15.5) in a form that are the same 
equations as those used to determine the same parameters of the C = C -bond MO of 
ethylene (Eqs. (14.242-14.268)) while matching the boundary conditions of the structure 
of benzene. 

164. The system of claim 161, wherein hybridization with 25% electron donation to each 
C = C -bond gives rise to the C bemene 2sp 3 HO-shell Coulombic energy 



3**- (CL-AP 1 ) given by Eq. (14.245). 



1 65. The system of claim 161, wherein to meet the equipotential condition of the union of 



3e 



the six C2sp 3 HOs, c 2 and C 2 of Eq. (15.42) for the aromatic C=C -bond MO is given 
by Eq. (15.62) as the ratio of 15.95955 eV , the magnitude of E QMb (c^^Zsp 3 ) (Eq. 

(14.245)), and 1 3.605804 eV , the magnitude of the Coulombic energy between the 
electron and proton of H (Eq. (1.243)): 

/ \ , , \ 13.605804 eV A 

C 2 (benzeneC2sp 3 HO)= c 2 (benzeneC2sp 3 HO)= - 0.85252 (15.143) 

wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to + 10%, if desired. 
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166. The system of claim 161, wherein the energies of each C=C bond of benzene are 
determined using the same equations as those of ethylene (Eqs. (14.251-14.253) and 
(14.319-14.333) with the parameters of benzene. 

u 

167. The system of claim 161, wherein ethylene serves as a basis element for the C=C 

bonding of benzene wherein each of the six C= C bonds of benzene comprises 
(0.75)(4)= 3 electrons according to Eq. (15.142). 

168. The system of claim 161, wherein the total energy of the bonds of the eighteen 
electrons of the C= C bonds of benzene, E T (q 6 H^C=c\ is given by (6)(o.75> times 

3e 

E T+Qse (C = C) (Eq. (14.492)), the total energy of the C=C -bond MO of benzene 
including the Doppler term, minus eighteen times E(C,2sp 3 ) (Eq. (14.146)), the initial 

energy of each C2sp* HO of each C that forms the C=C bonds of bond order two. 

169. The system of claim 161, wherein the total energy of the.six C=C bonds of benzene 
with three electron per aromatic bond given by Eq. (14.493) is 

■ 

E 7 [c 6 H 6 ,C=cj = (6X0.75^ (c = C>- (6X3>(c,2*p 3 ) 

= (6X0.75X"66.05796 eV}- 18(-14.63489 eV) (15.144) 
= -297.26081 eV - (-263.42798 eV) 
= -33.83284 eV 

wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to + 10%, if desired. 

1 70. The system of claim 161, wherein the results of benzene can be generalized to the 
class of aromatic and heterocyclic compounds. E hv of an aromatic bond is given by 
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E T (H 2 ) (Eqs. (1 1.212) and (14.486)), the maximum total energy of each # 2 -type MO 



such that 



S« = *fe> + *«)="k -31.63536831 eV 



\ 




(15.145), 



wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to ± 10%, if desired. 

171 . The system of claim 170, wherein the factor of 0.75 corresponding to the three 
electrons per aromatic bond of bond order two given in the Benzene Molecule ( C 6 H 6 ) 

r 

section modifies Eqs. (15.52-15.56). 



172. The system of claim 171, wherein the multiplication of the total energy given by Eq. 
(15.55) by f x as 0.75 with the substitution of Eq. (15.145) gives the total energy of the 
aromatic bond: 



wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to ± 10%, if desired. 

173. The system of claim 161, wherein the total bond energy of the aromatic group 




(15.146) 
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E D (group)— — 



f E(basis energies) + E T (atom - atom, msp 3 AO) 



fx 



W 



-31.63536831 eV 




Ane R 3 



m 



m t c 



3 r 3 



I" (f A« ( A0 1 c sKm (fs AO I HO)) 



) 



(15.147), wherein the calculated and measured values and constants recited in the 
equations herein can be adjusted, for example, up to + 10%, if desired. 

174. The system of claim 161, wherein benzene is considered as comprising chemical 
bonds between six CH radicals wherein each radical comprises a chemical bond between 
carbon and hydrogen atoms, energy components of V B , V p , T 9 V m , and E T are the same 

as those of the hydrogen carbide radical, except that E T (c = C,2sp 3 }= -1.13379 eV (Eq. 
(14.247)) is subtracted from E T (CH) of Eq. (13.495) to match the energy of each 

C- 77 -bond MO to the decrease in the energy of the corresponding C2sp 3 HO, wherein 

the calculated and measured values and constants recited in the equations herein can be 
adjusted, for example, up to ± 10%, if desired. 

175. The system of claim 174, wherein in the corresponding generalization of the aromatic 
CH group, the geometrical parameters are determined using Eq. (15.42) and Eqs. (15.1- 

15.5) with E T (atom - atom,msp l \AO)= -1 .13379 eV wherein the calculated and 

measured values and constants recited in the equations herein can be adjusted, for 
example, up to + 10%, if desired. 



176. The system of claim 174, wherein the total energy of the benzene C - H -bond MO, 
E Tt __ [C-H)> given by Eq. (14.467) is the sum of 0.5£ 7 (c = C,2sp* ), the energy 



change of each C2sp 3 shell per single bond due to the decrease in radius with the 



3c 



formation of the corresponding C= C -bond MO (Eq. (14.247)), and E T (ch\ the a 

benzene \ J 

MO contribution given by Eq. (14.441). In the corresponding generalization of the 
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aromatic CH group, the energy parameters are determined using Eqs. (1 5.146-15.147) 

(\ -1 13379 eV 
atom^atom 9 msp*.AO)= — : , wherein the calculated and 

measured values and constants recited in the equations herein can be adjusted, for 
example, up to + 10%, if desired 

177. The system of claim 1 76, wherein the energy contribution to the single aromatic CH 

3e 

bond is one half that of the C= C double bond contribution, which matches the energies 

3* 

of the CH and C-C aromatic groups, conserves the electron number with the equivalent 
charge density as that of 5 - 1 in Eqs. (15.18-15.21), and further gives a minimum energy 
for the molecule. 



3e 

178. The system of claim 177, wherein the breakage of the aromatic C-C bonds to give 
CH groups creates unpaired electrons in these fragments that corresponds to c = 1 in Eq. 



(15.56) with E mag given by Eq. (15.58). 



1 79. The system of claim 178, wherein each of the C - H bonds of benzene comprises two 
electrons according to Eq. (14.439). 

180. The system of claim 179, wherein from the energy of each C - H bond, 

-B D ( 12 C# ) (Eq. (14.477)), the total energy of the twelve electrons of the six C - H 

bonds of benzene, E r (C 6 H 69 C - H), given by Eq. (14.494) is 

E T (C 6 H 69 C- Hy (S^E^ ( l2 C//))= 6 (-3.90454 eV)= -23.42724 eV 

(15.148), wherein the calculated and measured values and constants recited in the 
equations herein can be adjusted, for example, up to + 10%, if desired. 

181. The system of claim 1 74, wherein the total bond dissociation energy of benzene, 

( u y 

E D (C 6 H 6 y given by Eq. (14.495) is the negative sum of E T \C 6 H 6 ,C=Cj (Eq. 
(14.493)) and E T (C 6 H S ,C - H) (Eq. (14.494)): 
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= -((-33.83284 eV)¥ (-23.42724 evty (15.149) 
= 57.2601 eV 

wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to + 10%, if desired. 

1 82. The system of claim 181, wherein using the parameters given in Tables 1 5 .2 1 4 and 
15.216 in the general equations (Eqs. (15.42), (15.1-15.5), and (15.146-15.147)) 
reproduces the results for benzene given in the Benzene Molecule (C 6 H 6 ) section as 

shown in Tables 15.214 and 15.216. 

» • 

183. A system of computing and rendering the nature of at least one specie selected from a 
group of diatomic molecules having at least one atom that is other than hydrogen, 
polyatomic molecules, molecular ions, polyatomic molecular ions, or molecular radicals, 
or any functional group therein, comprising physical, Maxwellian solutions of charge, 
mass, and current density functions of said specie, said system comprising: 

processing means for processing physical, Maxwellian equations representing charge, 
mass, and current density functions of said specie to produce at least one Maxwellian 
solution; and 

an output means for outputting the Maxwellian solution.' 

1 84. The system of claim 1 83, further comprising a data transfer system for inputting 
numerical data into or out of a computational components and storage components of the 
main system. 

185. The system of claim 183, further comprising a spreadsheet containing solutions of the 
bond parameters with output in a standard spreadsheet format. 

1 86. The system of claim 1 85, further comprising a data-handling program to transfer data 
from the spreadsheets into the main program. 
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1 87. The system of claim 1 83 9 wherein output may be at least one of graphical, simulation, 
text, and numerical data. 

188. The system of claim 1 87, wherein the output may be the calculation of at least one of: 
(1) a bond distance between two of the atoms; (2) a bond angle between three of the 
atoms; (3) a bond energy between two of the atoms; (4) orbital intercept distances and 
angles; and (5) charge-density functions of atomic, hybridized, and molecular orbitals, 
wherein the bond distance, bond angle, and bond energy are calculated from physical 
solutions of the charge, mass, and current density functions of atoms and atomic ions, 
which solutions are derived from Maxwell's equations using a constraint that a bound 
electron(s) does not radiate under acceleration 

1 89. The system of claim 1 83, wherein the charge, current, energy, and geometrical 
parameters are output to be inputs to other programs that can be used in further 
applications. 

190. The system of claim 189, wherein the data of heats of formation can be input to 
another program to be used to predict stability (existence of compounds) equilibrium 
constants and to predict synthetic pathways. 

m 

191 . The system of claim 190, wherein novel composition of matters can be discovered by 
calculating at least one of a bond distance between two of the atoms, a bond angle 
between three of the atoms, and a bond energy between two of the atoms, orbital intercept 
distances and angles, charge-density functions of atomic, hybridized, and molecular 
orbitals, the bond distance, bond angle, and bond energy being calculated from physical 
solutions of the charge, mass, and current density functions of atoms and atomic ions, 
which solutions are derived from Maxwell's equations using a constraint that a bound 
electron(s) does not radiate under acceleration. 

1 92. The system of claim 1 90, wherein the charge and current density functions can be 
used to predict the electric and magnetic fields of the species to determine other properties 
due to the interaction of the fields between species. 
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1 93 . The system of claim 1 92, wherein finite-element analysis is used to predict or 
calculate the interaction and resulting properties, such as the freezing point, boiling point, 
density, viscosity, and refractive index. 

194. The system of claim 1 83, wherein the output data can be used to give thermodynamic, 
spectroscopic, and other properties, aid in drug design and other applications with or 
without direct visualization. 

1 95 . The system of claim 1 83, wherein the data can be input into other programs of the 
system, which calculate thermodynamic and other properties, or performs a simulation, 
such as a chemical reaction or molecular dynamics. 

■ 

■ * 

1 96. A method of using any of the systems of claims 1 - 1 95. 

197. A use of any system of claims 1-195. 

1 98. A composition of matter comprising a plurality of atoms, the improvement comprising 
a novel property or use discovered by calculation of at least one of 

a bond distance between two of the atoms, 
a bond angle between three of the atoms, and 
a bond energy between two of the atoms, 
orbital intercept distances and angles, 

charge-density functions of atomic, hybridized, and molecular orbitals, 
the bond distance, bond angle, and bond energy being calculated from physical solutions of 
the charge, mass, and current density functions of atoms and atomic ions, which solutions are 
derived from Maxwell's equations using a constraint that a bound electron(s) does not radiate 
under acceleration* 

1 99. A composition of matter according to claim 1 98, wherein the novel property is a new 
pharmaceutical use. 

200. A composition of matter according to claim 1 98 that is novel, wherein the novel 
property is stability at room temperature of a new arrangement of atoms or ions. 
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201 . A novel composition of matter discovered by calculating at least one of 
a bond distance between two of the atoms, 

a bond angle between three of the atoms, and 
a bond energy between two of the atoms, 
orbital intercept distances and angles, 

charge-density functions of atomic, hybridized, and molecular orbitals, 
the bond distance, bond angle, and bond energy being calculated from physical solutions of 
the charge, mass, and current density functions of atoms and atomic ions, which solutions are 
derived from Maxwell's equations using a constraint that a bound electron(s) does not radiate 
under acceleration* 

■ ... 

202. A system of determining at least one of 
a bond distance between two of the atoms, 

a bond angle between three of the atoms, and 
a bond energy between two of the atoms, 
orbital intercept distances and angles, 

charge-density functions of atomic, hybridized, and molecular orbitals, 
wherein the bond distance, bond angle, and bond energy being calculated from physical 
solutions of the charge, mass, and current density functions of atoms and atomic ions, which 
solutions are derived from Maxwell's equations using a constraint that the bound electron(s) 
does not radiate under acceleration. 

203. The system of claim 202, further comprising discovering a new composition of 
matter. 

> 

204. A system of claim 1-195, wherein the calculated and measured values and constants 
recited therein can be adjusted up to + 1 0% as desired. 
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Fig. 3 
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